f"««*EEE 


\ 


MCROOOnr  RESOLUTION  TEST  CHART 
MTHNN.  MJAEAU  Of  STANOAMS-UO-A 


>4  REPORT  C UNITY  CLASSIFKA 

UNCLASSIFIED 


eSECURiTV  classification  authority 


2b  OECLASSiFiCATiON  / DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


REPORT  DOCUMENTATION  PAGE 


lb  RESTRICTIVE  MARKINGS 


)  DISTRIBUTION/ AVAILABILITY  0f  REPORT 

Approved  for  public  release; 
Distribution  is  unlimited 


S  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


64  NAME  OF  PERFORMING  ORGANIZATION  6b  OFFICE  SYMBOL  7e  NAME  OF  MONITORING  ORGANIZATION 

(if  ipplkibl*) 

Naval  Postgraduate  School  55  Naval  Postgraduate  School 


6c  ADORE  SS  (City.  Stit*.  end  IIP  Cod*) 

Monterey,  CA  93943-5000 


7b  ADDRESS  (City.  Stit*.  and  it*  Cod*) 

Monterey,  CA  93943-5000 


Ba  NAME  OF  FUNOING/ SPONSORING 
ORGANIZATION 


8b  OFFICE  SYMBOL 
(If  ipplkibl*) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


Be  ADDRESS  (City.  Stitt.  *nd  IIP  Cod*) 


ii  Title  ( Include  Security  Cliuifkition) 

THE  CURVE  FITTING  OF  PORTABLE  WEAR  MJ 


iz  PERSONAL  AUTHOR(S)  X 

MIN,  Byuna  Ho  / 


'6  SUPPLEMENTARY  NOTATION 


COSATi  COOES 


GROUP  I  SUBGROUP 


1  10  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM 
ELEMENT  NO 

PROJECT 

NO 

TASK 

NO 

WORK  UNIT 
ACCESSION  NO 

LL  ANALYSIS 


14  DATE  OF  REPORT  (Y**r.  Month.  Oiy)  IS  PAGE  COoXl 

1987  March  67  X 


iByUBJECT  TERMS  (Continue  pn  revert*  if  nutuiry  *nd  identify  by  block  number) 


Curvei.Fi  t  ting  j 
Equation^  spectrosc 


,  Rational,*. 


copy . 


TRACT  (Continue  on  reverie  if  ntttuiry  ind  i dtntify  by  block  number) 

This  thesis  describes  the  curve  fitting  algorithm  that  the  Portable 
Wear  Metal  Analyzer  used  for  calculating  concentrations  in  ppm  and 
compares  this  with  some  alternative  algorithms.  Each  algorithm 
considered  fits  a  curve  to  the  three  standards  used,  which  were  20%, 
50%,  and  100%  of  the  full  scale  for  all  nine  wear  metals  APL  was 
used  for  all  programs.  xs  /)  " 


;0  D  S*R'3UTiON/ AVAILABILITY  of  ABSTRACT 
fDuNCLASSlFiEO'UNLiMlTEO  □  SAME  AS  RPT 

□  OTIC  USERS 

Hi  name  OF  RESPONSIBLE  iNOiVIOUAL 

21  ABSTRACT  SECURITY  CLASSIFICATION 


OO  FORM  1473, 84  mar 


8)  APR  edition  miy  b *  vied  until  e«h»u»ted 
All  other  edition!  ere  obsolete 


SECURITY  CLASSIFICATION  of  This  page 

UNCLASSIFIED 


Approved  for  public  release;  distribution  is  unlimited. 


;.v  The;  Curve  Fitting 
of 

Portable  Wear  Metal  Analysis 


by 


Byung  II.  Min 
Major,  Korean  Army 
B.S.,  Korea  Militarv  Academy,  1976 
B.S.,  Seoul  National  University,  1980 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
March  1987 


2 


ABSTRACT 


This  thesis  describes  the  curve  fitting  algorithm  that  the  Portable  Wear  Metal 
Analyzer  used  for  calculating  concentrations  in  ppm  and  compares  this  with  some 
alternative  algorithms.  Each  algorithm  considered  fits  a  curve  to  the  three  standards 
used,  which  were  20%,  50%,  and  100%  of  the  full  scale  for  all  nine  wear  metals.  APL 
was  used  for  all  programs. 
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I.  INTRODUCTION 


Spectrometric  oil  analysis  is  used  to  determine  the  type  and  amount  of  wear 
metals  in  lubricating  oil  samples.  The  Portable  Wear  Metal  Analyzer  (PWMA)  was 
developed  for  the  USAF  to  perform  engine  oil  analysis  when  aircraft  are  deployed 
away  from  their  home  bases.  The  analyzer  can  detect  nine  wear  metals  in  parts  per 
million  concentration  ranges.  PWMA  takes  samples  with  known  concentration  to 
obtain  absorbance  areas,  and  uses  those  to  determine  a  calibration  curve.  Absorbance 
areas  obtained  from  unknown  samples  are  compared  to  the  calibration  curve  so  that 
their  concentration  can  be  calculated. 

This  thesis  describes  the  curve  fitting  algorithm  that  the  PWMA  used  for 
calculating  the  concentration  in  ppm  and  compares  this  with  some  alternative 
algorithms  for  the  curve  fitting.  We  compare  these  algorithms  to  see  which  algorithm 
performs  best.  Each  algorithm  considered  fits  a  curve  to  the  three  standards  used  20%, 
50%,  and  100%  of  the  full  scale  for  all  nine  wear  metals.  APL  was  used  for  all 
programs. 

Data  gathered  as  the  result  of  the  field  test  of  the  PWMA  are  described  in 
Chapter  II.  Three  rational  functions,  and  polynomials  which  are  quadratic  in 
absorbance  and  quadratic  in  concentration  are  suggested  as  alternatives  in  Chapter  III. 
A  squared  error  technique  is  employed  to  compare  the  algorithms'  performances. 


II.  PWMA  AND  DATA  DESCRIPTION 

A.  GENERAL 

The  PWMA  test  program  was  conducted  at  Naval  Air  Station,  Pensacola  FL, 
and  three  Air  Force  field  test  sites,  Langley  AFB,  VA,  Myrtle  Beach  AFB,  SC,  and 
Elmendorf  AFB,  AK  in  1985.  The  purpose  of  this  test  plan  was  to  determine  the 
effectiveness  and  evaluate  the  functionality  of  the  newly  developed  PWMA  prototypes 
for  measuring  wear  metal  content  in  turbine  engine  lubricating  oils  under  field 
conditions. 

When  aircraft  are  deployed  from  their  home  bases  to  austere  sites  where  normal 
oil  analysis  capabilities  are  nonexistent,  there  still  exists  a  requirement  to  determine 
wear  metal  concentration  in  turbine  engine  lubricants.  Due  to  the  critical  nature  of 
turbine  engine  lubrication  system  component  wear,  particularly  in  high  performance 
systems  such  as  the  F-15,  F-16,  and  A-10  aircraft,  the  requirement  for  oil  analysis  may 
be  on  an  "after  each  flight"  basis.  Presently,  the  Air  Force  deploys  Flame  Atomic 
Absorption  (FAA)  spectrophotometers  [Ref.  1:  p.  50],  which  were  not  designed  for 
mobile  use.  These  instruments  are  sensitive  and  highly  susceptible  to  damage  when 
deployed  for  mobility  exercises;  such  damage  may  lead  to  partial  or  complete  loss  of  oil 
analysis  capability  in  early  stages  of  such  exercises.  In  addition,  the  instruments  require 
hazardous  support  accessories,  such  as  compressed  gases  (nitrous  oxide  and  acetylene) 
and  solvents  for  dilution  (and,  in  some  past  cases,  acid  solutions  for  digestion  of  wear 
metal  particulates). 

In  1980  FIQ  TAC  issued  TAF  SON  305-80,  establishing  the  requirement  for  a 
critically  needed  light-weight  portable  aircraft  engine  oil  analysis  unit.  In  response, 
AFWAL1  undertook  a  contractual  effort  to  determine  the  feasibility  of  developing  such 
a  device  (Contract  Number  F33615-80-C-2037).  Since  no  such  spectrometer  with 
required  light-weight  portable  features  was  commercially  available,  a  major 
development  program  was  initiated  to  ruggedize  and  miniaturize  the  conceptual 
hardware  for  mobility  use.  The  actualization  of  the  portable  wear  metal  analyzer  has 
been  achieved  under  Contract  Number  F33615-81-C-2080  by  Perkin-Elmer 
Corporation,  Applied  Science  Division. 


^ir  Force  Wright  Aeronautical  Laboratory 
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B.  THE  PORTABLE  WEAR  METAL  ANALYZER 

The  PWMA  features  simultaneous  measurement  of  nine  wear  metals  in 
lubricating  oils  with  an  analysis  cycle  of  four  minutes  (via  air  cooling):  Silver  (Ag), 
Aluminum  (Al),  Chromium  (Cr),  Copper  (Cu),  Iron  (Fe),  Magnesium  (Mg),  Nickel 
(Ni),  Silicon  (Si),  and  Titanium  (Ti).  Moreover,  with  the  potential  modification  to 
water  cooling  on  production  units,  this  time  can  be  reduced  to  two  minutes.  In 
addition,  the  PWMA  may  be  operated  on  either  110/220  VAC,  50  or  60  Hertz,  to 
accommodate  use  in  both  the  U.S  and  Europe.  The  PWMA  is  packed  in  two 
aluminum  “suitcases”:  the  electronics/argon  case,  weighing  forty  pounds;  and  the 
furnace;  optics  case,  weighing  sixty  pounds.  The  cases  are  readily  set-up  and 
interconnected  by  two  multiconductor  cables  and  one  pneumatic  hose.  Developed  for 
simple  operation  and  reduced  dependence  on  skill  and  labor  of  the  operator,  the 
PWMA  is  characterized  as  follows:  ease  of  set-up  and  self-contained  configuration; 
minimal  consumables  (sample  tips,  graphite  tubes  and  argon  gas  supply);  no  hazardous 
gases  or  chemicals  (inert  compressed  gas  only);  microprocessor-controlled  software  for 
calibration,  diagnostic  etc.;  and  elimination  of  transcription  errors  by  automated  read¬ 
out/print-out  of  data  and  results. 

The  PWMA  operates  on  the  principle  of  absorption  of  characteristic  radiation 
by  atoms  of  the  wear  metals  of  interest.  Two  multi-element  hollow  cathode  lamps 
(HCLs),  alternately  pulsed  for  source  combination  and  modulation,  provide  emission 
source  lines  for  discrete  line  absorption  by  specific  atoms,  which  are  formed  through 
heating  in  a  preprogrammed  ramped  cycle  (up  to  3000°  C)  in  the  graphite  furnace.  All 
nine  absorption  lines  can  be  separated  simultaneously  by  the  polychrometer  and 
transmitted  to  respective  discrete  photomultiplier  tubes  (where  the  signal  is  transduced 
from  energy  to  current).  Next,  each  individual  electrometer  converts  the  current  to  a 
voltage  signal,  which  is  amplified  and  processed  by  the  analog  (A/D)  converter.  The 
output  is  printed  in  concentration  (parts  per  million,  ppm)  for  all  measurements. 

The  operation  of  the  PWMA  is  as  follows.  The  operator  first  introduces  an  oil 
sample  into  the  graphite  tube  of  the  Graphite  Furnace  Assembly.  The  PWMA  start 
button,  located  on  the  control  panel,  is  then  depressed,  and  the  entire  analysis  program 
is  initiated  and  controlled  by  an  "on  board"  microcomputer.  As  the  graphite  furnace 
executes  the  atomization  cycles,  the  emission  line  spectra  originating  from  the  HCLs 
are  focused  through  the  graphite  tube  to  the  polychrometer  where  they  are  separated 
into  individual  spectral  lines  and  are  further  conveyed  to  the  exit  slits.  The  spectral 
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lines  pass  through  their  respective  exit  slits  where  their  spectral  intensities  are  detected 
and  measured  by  individual  Photomultiplier  Tubes  (PMTs).  Wear  metals  present  in 
the  oil  sample  will  absorb  the  light  associated  with  each  metal’s  spectral  line,  thus 
reducing  the  intensities  reaching  the  PMTs.  The  amount  of  light  absorbed  is  related  to 
the  concentration  of  metal  in  the  oil  sample. 

The  PWMA  has  several  function  modes.  In  the  calibration  mode,  a  prepared 
sample  with  known  concentrations  is  used  to  get  the  absorbance  data.  Analysis  mode  is 
selected  to  analyze  the  unknown  samples.  The  prototype  PWMA  prints  out  two  types 
of  data  in  the  analysis  mode,  concentration  data  (computed  by  the  on-board 
microprocessor)  and  absorbance  data  (in  parentheses).  These  data  can  be  used  to 
compare  analysis  algorithms. 

C.  TEST  PROGRAM  AND  DATA  COLLECTED 

The  four  prototypes  were  tested  at  their  respective  field  evaluation  sites.  The  test 
followed  the  PWMA  field  test  bum  sequence,  which  is  in  Appendix  A.  The  graphite 
tubes,  used  in  the  graphite  furnace,  deteriorate  gradually  with  use,  and  should  be 
replaced  after  every  160  sample  bums.  If  a  new  graphite  tube  is  being  used, 
preconditioning  of  the  tube  is  required.  Laboratories  were  to  follow  the  test  sequence 
which  included  calibration  samples,  verification  samples,  correlation  samples,  random 
samples,  reslope,  and  so  on.  Three  calibration  samples  were  prepared  in  MIX-L  7S0S 
oil  containing  20%,  50%  and  100%  of  full  scale  for  all  nine  metals.  Perkin-Elmer 
provided  these  samples.  Three  verification  samples  were  prepared  in  MIX-L  7808  oil 
containing  10%,  40%  and  70%  of  full  scale  for  all  nine  metals.  Perkin-Elmer  provided 
these  samples,  too.  Correlation  Samples  used  turbine  engine  oil  samples.  Each  sample 
contained  a  minimum  of  7  wear  metals.  Random  samples  were  routine  daily  samples 
obtained  at  each  evaluation  site.  This  analysis  employs  only  the  calibration  and 
verification  samples. 

The  concentration  data  printed  out  by  the  PWMA  are  calculated  by  the 
microprocessor  with  the  use  of  a  calibration  technique.  Samples  with  known 
concentrations  are  analyzed  first  and  the  absorbance  data  obtained  are  used  to 
construct  calibration  curves.  Absorbance  data  obtained  from  unknown  samples  are 
then  compared  to  the  calibration  curves  so  that  their  concentrations  can  be  calculated. 
Later  verification  samples  obtained  from  known  samples  give  an  indication  of  how 
valid  the  calibration  curve  is. 
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Table  1  shows  data  collected  from  one  of  the  field  test  sites.  P523  means 
"Pensacola  May  23rd,  1985",  the  place  and  the  date  the  data  were  collected,  "401"  is 
the  tube  sequential  number.  CAL  I A  and  CAL  IB  stands  for  the  calibration  sample 
20%.  Cal  II  and  Cal  III  are  50%  and  100%  of  full  scales  respectively.  In  order  to 
provide  a  more  accurate  result,  a  two-point  average  calibration  was  used  for  the  field 
test,  i.e.,  two  CAL  I  samples  were  run  at  CAL  I  position;  two  CAL  II  samples  at  CAL 
II  position;  and  two  CAL  III  samples  at  CAL  III  position.  For  each  calibration  run, 
the  numerical  values  are  the  absorbance  data  printed  out  (rounded  to  integers).  The 
PWMA  averaged  the  two  samples  at  the  same  level,  and  stored  them  as  the  calibration 
points,  but  the  average  value  was  not  printed  out. 

TABLE  1 

ABSORBANCE  DATA  OF  TUBE  401 


P523  401 

Element 

Fe 

Ag 

A1 

Cr 

Cu 

Mg 

Ni 

Si 

Ti 

CAL  IA 

58 

53 

10 

38 

23 

44 

40 

10 

10 

CAL  IB 

56 

52 

10 

34 

22 

46 

38 

8 

9 

CAL  IIA 

115 

98 

19 

83 

47 

60 

92 

28 

30 

CAL  I IB 

111 

100 

21 

84 

43 

61 

96 

18 

29 

CAL  I IIA 

143 

129 

26 

104 

70 

72 

119 

33 

42 

CAL  I I IB 

145 

131 

24 

104 

71 

74 

120 

35 

44 

VERIF  10 

CONC 

9 

1 

2 

1 

4 

3 

3 

2 

5 

VERIF  10 

ABS 

28 

28 

5 

17 

12 

32 

19 

4 

11 

VERIF  40 

CONC 

44 

4 

10 

4 

18 

9 

12 

9 

7 

VERIF  40 

ABS 

105 

90 

17 

70 

42 

56 

79 

20 

20 

VERIF  70 

CONC 

70 

7 

9 

6 

28 

17 

20 

12 

10 

VERIF  70 

ABS 

129 

114 

16 

88 

56 

66 

105 

26 

29 

Part  of  the  data  for  tube  401  is  shown  in  Table  1.  The  complete  data  for  tube 
401  is  attached  in  Appendix  B.  VERIF  10,  VERIF  40  and  VERIF  70  are  the 
verification  samples.  The  data  in  the  first  row  of  VERIF  10  is  the  computed 
concentration  and  the  second  row  is  absorbance  data. 


Table  2  shows  all  the  data  available  for  this  study,  Langley  by  location  tested  6 
tubes.  Myrtle  Beach  tested  5  tubes,  Pensacola  tested  1 1  tubes,  and  Elmendorf  tested  6 
tubes.  All  28  tubes  were  used  during  the  field  test  program. 

TABLE  2 

DATA  OF  EACH  LABORATORY 


Lab. 

Tubes 

Total 

Langley 

FPM201 

PPM204 

M 

PPM203 

PPM207 

6 

Myrtle 

Beach 

PPM302 

PPM305 

PPM303 

PFM306 

FPM304 

5 

Pensacola 

PPM401 

PPM405 

PFM409 

PFM415 

PPM403 

PPM406 

PPM412 

FPM416 

PPM404 

PPM408 

PPM413 

11 

Elmendorf 

PPM501 

PPM504 

PPM502 

PFM505 

PPM503 

PPM507 

6 

In  the  PWMA,  three  calibration  standards  of  20%,  50%  and  100%  are  used  to 
construct  the  calibration  curve.  The  absorbance  data  in  CAL  IA  through  CAL  IIIB 
are  used  to  compute  the  calibration  curve  coefficients.  A  pair  of  two  point  rational 
equations  were  used  to  fit  the  calibration  curves,  one  for  each  of  the  nine  elements. 
The  function  used  is  given  in  equation  2.1. 

C  -  K,A/(K2A-1)  (2.1) 

where  Kj,  K2  are  unknown  coefficients,  C  is  concentration  in  ppm,  and  A  is 
absorbance  data. 

The  first  part  of  the  curve  uses  absorbance  values  obtained  from  20%  and  50% 
samples  and  the  second  part  of  the  curve  uses  values  obtained  from  50%  and  100%. 
We  will  explain  how  the  concentration  value  of  9,  the  first  line  of  the  VERIF  10  in 
Table  1,  was  computed  in  the  PWMA.  Figure  2.1  shows  the  calibration  technique 
used  by  the  PWMA.  In  Figure  2.1,  the  first  part  of  the  curve  for  the  FE  is  given  by 
the  solid  line,  the  second  part  of  the  curve  is  given  by  the  dashed  line,  and  symbol  + 
identifies  the  three  calibration  points.  The  verification  absorption  values  of  129,  105 


Figure  2.1  Calibration  Technique  used  by  the  PWMA. 

and  28  are  applied  to  the  curve,  and  the  calibrated  data  computed  by  the  PWMA  are 
70,  44  and  9  respectively.  However,  the  correct  values  are  70,  40  and  10.  Verification 
sample  VERIF  10  which  has  known  concentration  10  produces  absorbance  data  28,  as 
shown  in  Table  1.  If  the  absorbance  data  28  was  fitted  to  the  curve  from  the 
absorbance  axis,  then  its  concentration  should  be  close  to  10.  However,  it  show's  some 
deviation  because  the  curve  point  projected  to  the  concentration  axis  makes 
concentration  9,  as  shown  in  Figure  2.1.  When  the  observed  absorbance  is  within  the 
range  of  0  to  the  absorbance  value  at  50%  concentration,  the  first  curve  is  applied  to 
compute  the  concentration.  For  the  absorbance  data  greater  than  this,  the  second 
curve  is  applied  to  compute  the  concentration.  The  computer  selects  the  proper  curve 
to  calculate  data,  dependent  upon  the  calibrating  concentrations  of  all  nine  elements. 
(Ref.  2:  p.  67) 

The  PWMA  has  a  multi-element  capability,  the  9  elements  already  listed  in  Table 
1,  Iron,  Silver,  Aluminum,  Chromium,  Copper,  Magnesium,  Nickel,  Silicon  and 
Titanium. 
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The  full  scale  capabilities  for  these  elements  are  shown  in  Table  3.  For  example, 
the  full  scale  value  for  Fe  is  100,  and  is  10  for  Ag.  The  columns  labeled  20%,  50%  and 
100%  are  CAL  I,  CAL  II  and  CAL  III  samples  respectively.  The  10%,  40%  and  70% 
columns  are  VERIF  10,  VERIF  40  and  VERIF  70  respectively.  The  concentration 
varies  with  the  sample  and  the  element. 

TABLE  3 

TABLE  OF  OIL  SAMPLE  MIXTURES 


METAL 

concentration  (ppm) 

10% 

20% 

40% 

50% 

70% 

100% 

FE 

10 

20 

40 

50 

70 

100 

AG 

1 

2 

4 

5 

7 

10 

AL 

2.  5 

5 

10 

12.  5 

17.5 

25 

CR 

1 

2 

4 

5 

7 

10 

CU 

4 

8 

16 

20 

28 

40 

MG 

2.5 

5 

10 

12.  5 

17.5 

25 

NI 

3 

6 

12 

15 

21 

30 

SI 

2 

4 

8 

10 

14 

20 

TI 

2 

4 

8 

10 

14 

20 

CAL  I 

CAL  II 

CAL  III 

III.  ALGORITHMS 


A.  GENERAL 

If  <me  increment  of  concentration  C  causes  an  equal  increment  of  absorbance  A 
at  all  levels  of  C,  then  we  say  that  there  is  no  curvature  in  the  analytical  curve 
expressing  the  relation  between  C  and  A.  In  reality,  in  the  presence  of  stray  light  the 
analytical  curve  bends  towards  the  concentration  axis.  [Ref.  3:  p.  681] 

Calibration  curvature  can  be  controlled  to  an  extent  by  the  instrument  design. 
The  degree  of  curvature  for  iron  at  the  248.3-nm  wavelength,  for  example,  is  partially 
controlled  by  the  spectral  bandwidth  since  most  of  the  curvature  can  be  attributed  to 
overlap  of  nearby  nonabsorbed  lines  in  the  iron  spectrum.  This  curvature  effect  is 
analogous  to  stray  light  interferences  observed  in  molecular  spectroscopy.  However,  at 
practical  instrument  settings,  it  is  not  possible  to  eliminate  this  curvature  completely. 

B.  E.  Limbeck  investigated  transformations  of  the  data  from  AA  (Atomic 
Absorption)  instruments  to  forms  that  would  produce  more  easily  handled  curves. 
Graphic  trials  on  several  possible  transformations  revealed  that  a  plot  of  A/C  vs  A  or 
C/A  vs  C  produced  pseudocalibrations  that  were  much  less  curved  than  the  A  vs  C 
calibrations,  even  for  the  extreme  cases.  He  concluded  that  over  a  limited 
concentration  range,  two  or  three  standards  were  sufficient  to  produce  excellent  fits  for 
AA  calibrations  with  the  rational  algorithm.  [Ref.  4:  p.  96] 

As  already  mentioned,  the  PWMA  uses  a  rational  function  which  uses  two 
standards;  to  cover  the  full  range  this  equation  is  used  twice.  This  same  function  will 
be  used  twice  in  the  following  analysis.  The  PWMA  on  board  computer  uses  single 
precision  arithmetic  and  establishes  its  calibration  curve  internally.  The  data  available 
for  this  study  was  printed  out  in  the  PWMA  paper  tapes,  which  rounded  the  observed 
absorbance  numbers.  Thus,  if  one  uses  the  PWMA  equation  to  establish  a  curve,  and 
then  reads  the  rounded  absorbance  figures  back  through  this  curve,  one  does  not 
exactly  arrive  at  the  concentration  printed  out.  Thus  the  first  algorithm  should  give 
essentially  the  same  results  as  those  labeled  PWMA,  which  were  the  rounded  figures 
from  the  internal  curve.  The  other  two  rational  functions  have  three  unknown 
coefficients  so  that  they  are  used  one  time  with  three  calibration  standards.  We  will 
examine  three  rational  functions  and  two  quadratic  functions,  which  are  defined  below. 


Algorithm  1  :  C  -  A  /  (  K,  +  K2  A  ) 


Algorithm  4 :  C 


Algorithm  5 :  A 


(  K,  A  +  K3A2)/(  K2A-  1) 
A  /  (  Kj  +K2A  +  K}  A2) 

K,  +  K2A  +  K3A2 
K1  +  K2  C  +  K3  C2 


B.  ALGORITHM  1 

This  is  useful  when  a  simple  linear  approximation  gives  a  good  fit  to  the  plot  of 
A/C  versus  A.  The  following  rational  function  was  chosen  as  algorithm  1: 

A/C-Kl  +  K2A  (3.1) 

It  is  clear  by  a  simple  rearrangement  of  the  formula  that  the  concentration  is  a 
function  of  absorbance: 


C  -  A  /  (  Kj  +  K2  A  ) 


(3.2) 


For  brevity,  the  calibration  algorithm  using  this  linear  rational  function  is 
referred  to  as  the  'rational  method'  and  is  the  one  actually  employed  in  the  PWMA. 
This  algorithm  uses  a  pair  of  two-point  rational  equations  3.2  to  fit  data  for  all  nine 
elements.  The  first  part  of  the  curve  uses  absorbance  data  obtained  from  20%  and 
50%  samples,  and  the  second  part  of  the  curve  uses  data  obtained  from  50%  and 
100%  samples. 

C.  ALGORITHM  2 

William  B.  Barnett  found  that  the  equation. 


C 


K0- 


K,  A  +  K3  A2 
K2  A  -  1 


(3.3) 


was  applicable  in  fitting  the  analytically  useful  concentration  range  of  atomic 
absorption  calibration  curves.  [Ref.  5:  p.  829J 
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In  tins  equation,  C  is  the  concentration,  A  is  the  observed  absorbance,  and  the 
Kj's  are  coefficients  which  must  be  determined  during  calibration.  The  coefficient  K0  is 
the  reslope  constant  which  is  set  equal  to  1.0  during  the  initial  standardization,  which 
gives 


C 


K,  A  +  K,  A2 
K2  A  -  1 


(3.4) 


In  equation  3.4,  concentration  has  been  expressed  as  a  function  of  absorbance.  This  is 
the  reverse  of  the  usual  representation  where  absorbance  is  considered  to  be  a  function 
of  the  standard  concentrations.  However,  expressing  concentration  as  the  dependent 
variable  facilitates  its  computation  without  the  need  for  extracting  roots  of  equations 
as  others  had  done.  Equation  3.4  can  actually  be  used  in  three  forms:  as  shown  in  the 
PWMA,  with  K3  being  set  equal  to  zero, 


C  -  K,  A  /  (  K2  A  -1  )  (3.5) 

or  with  K2  and  K3  both  being  set  equal  to  zero.  This  latter  case  is,  of  course,  the 
equation  for  a  straight  line  which  passes  through  the  origin. 

When  the  number  of  standards  matches  the  number  of  coefficients  to  be 
determined,  it  is  a  simple  matter  to  determine  the  coefficients  by  the  solution  of 
simultaneous  equations.  We  will  employ  simultaneous  equations,  to  show  how  to  get 
the  coefficients  in  APL. 

Suppose 

C  -  (  Cj,  c2,  c3  ),  A  •*  (  8j,  a2,  a3 ). 

If  we  express  both  sides  of  equation  3.4  in  terms  of  components  and  then  equate 
corresponding  components,  we  obtain  the  system 


c, 


Lh  ilL. 


K2a. 


1 


i  -  1,2,3 


(3.6) 


which  leads  to 


K2  Sj  ^ 


Ci  -  K1  ai  +K3ai2 


i  -  1,  2,  3 


(3.7) 
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i  -  I,  2,  3 


(3.8) 


Cj  -  -  Kj  »i  +  K2  aj  c,  -  K3  a.2 
Define 

K  -  (  Kr  K2,  Kj  ).  and  P  -  ||  py  || 

with 

Pi,  “  * 

Pi2* 

Pi3  “  *  ai2- 

Then 

C  -  P  K  (3.9) 

The  solution  for  K  is 
K  +  CMP 

where  I  is  the  matrix  inverse  operator  in  APL  (Ref.  6:  pp.  82,  83(. 

While  this  curve  may  well  faithfully  represent  the  Absorbance  -  Concentration 
relationship  over  the  useful  range,  there  are  certain  possible  observed  values  where  it 
does  not.  For  example  algorithm  2  will  not  perform  well  when  the  observed  pairs  (c., 
a.)  lead  to  K2  *  1/a,  where  at  <  a  <  aJf  since  the  denominator  of  equation  3.6 
becomes  zero.  For  example,  with  C  *  (  20,  50,  100  ),  A  **(  43,  117,  152  ),  the 
calibration  coefficients  K  are  (  -0.4693,  0.00694,  0.003322).  The  fitted  curve  defined  by 
these  coefficients  is  shown  in  Figure  3.1.  The  horizontal  line  is  an  asymptote  line  at  A 
equal  to  144.16,  which  is  the  reciprocal  of  K2-  If  the  observed  absorbance  is  slightly 
less  than  144.16,  the  computed  concentration  will  be  a  large  negative  number,  and 
when  the  absorbance  data  is  slightly  greater  than  144.16,  then  the  computed 
concentration  will  be  a  large  positive  number,  a  behavior  which  would  not  be 
acceptable. 
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Figure  3. 1  An  Asymptote  Line  in  Algorithm  2. 
D.  ALGORITHM  3 

This  algorithm  was  suggested  by  Limbeck  (Ref,  4:  p.  961, 


K,  +K2A  +  Kj  A* 


(3.10) 


This  is  a  dilfercnt  extension  of  algorithm  1,  namely,  A/C  •  K,  +  K2  A.  In 
equation  3.10,  the  quadratic  term  is  included.  To  get  the  coefficients  in  APL,  we  again 
employ  the  notation  of  simultaneous  equations. 


K,  +  K2a.  +  Kj  a.2  -  i  -  1.  2,  3 


(3.11) 


Now,  let 


C  « 
Pi! 


(  3,  'Cj,  a,  c2,  aj/Cj )  with 

l. 
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and  again  C  *  P  k,  so  the  solution  is 
K  +  CiP 

This  algorithm  can  also  produce  a  curve  with  an  interior  asymptote,  like 
algorithm  2,  when  the  denominator  is  zero. 

Another  type  of  undesirable  shape  can  occur  with  this  case  (as  well  as  algorithm 
2).  This  is  illustrated  by  taking 

C  -  (  20,  50,  100  ), 

A  -  (  43,  75,  152  ). 

Then,  we  find 

K  -  (  3.632,  -0.0425,  0.000188  ). 

The  resulting  curve  is  pictured  in  Figure  3.2.  Increased  absorbance  will  cause  lower 


Figure  3.2  Possible  Observed  Curve  from  Algorithm  3. 
estimated  concentration  for  A  >  140  with  this  example. 
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E.  ALGORITHM  4 

A  simple  quadratic  equation 

C  *  Kj  +  K2  A  +  K3  A2  (3.12) 

might  be  a  reasonable  alternative  for  a  calibration  curve. 

Equation  3.12  has  been  expressed  as  a  function  of  absorbance.  This  is  the  reverse 
of  the  usual  representation  where  absorbance  is  considered  to  be  a  function  of  the 
standard  concentrations.  However,  expressing  concentration  as  the  dependent  variable 
facilitates  its  calculation  without  the  need  for  extracting  roots  of  equations  as  is  needed 
for  algorithm  5. 

This  algorithm  may  give  a  satisfactory  shape,  depending  on  the  A  values  which 
occur  during  calibration  bums. 

It  has  an  undesirable  shape  for 

C  -  (20,  50,  100), 

A  -  (40,  130,  170), 

K  -  (43.3  -0.865,  0.007). 

Pictured  in  Figure  3.3,  here  decreasing  A,  for  A  <  60,  gives  increasing  C. 

When 

C  -  (20,  50,  100), 

A  -  (40,  60,  170), 

K  -  (-59.3,  2.3,  -0.008), 

algorithm  4  gives  the  curve  showii  in  Figure  3.4,  in  which  increasing  A,  for  A  >  140, 
leads  to  decreasing  C. 

Note  that  neither  of  these  two  curves  pass  through  the  origin,  since  the 
polynomial  equation  will  not  necessarily  give  zero  concentration  with  zero  absorbance. 
In  APL  the  coefficients  K  are  given  by 
K  +  CMA •  .  *  0 1 2 


figure  3.3  The  Curve  Bending  Downward  at  Lower  Concentration  by  Algorithm  4. 

where  C  -  (Cj,  c2,  )  are  the  concentrations  used  and  A  ■  (  aJt  a2,  a}  )  are  the 

observed  absorbance  numbers. 

F.  ALGORITHM  5 

The  final  curve  to  be  considered  is 

A  -  Kt  +K2C  +  KjC2  (3.13) 

In  equation  3.13,  absorbance  has  been  expressed  as  a  function  of  concentration 
which  is  the  natural  representation.  However,  expressing  absorbance  as  the  dependent 
variable  requires  the  solution  of  a  quadratic  equation  to  find  concentration.  Equation 
3. 1 3  can  be  expressed  like  this, 

K3C2  +  K2C  +  K,-A  -  0  (3.14) 

Solving  the  resulting  equation  for  C,  we  obtain  the  formula, 


Figure  3.4  The  Curve  Bending  Upward  at  High  Concentration  by  Algorithm  4. 


„  -K,±vrK,2-4K,(  K,  -  A  ) 

c - 2 - r...l - - —  (3.15) 

Whenever  K22  -  4  K3  (  K,  -  A  )  is  negative,  there  is  no  concentration  value  C  for  the 
given  absorbance  A. 

As  with  algorithm  4,  two  different  undesirable  shapes  may  occur;  if  the  initial 
calibration  data  is 

C  -  (20,  50,  100), 

A  -  (40,  140,  170), 

then 

K  -  (-60.8,  5.7,  -0.03), 

giving  the  curve  pictured  in  Figure  3.5. 

If  the  initial  data  is 


Figure  3.5  The  Curve  Bending  Downward  at  High  Concentration  by  Algorithm  5. 

C  -  (20,  50,100), 

A  -  (40,  50,  170), 

then 

K  -  (59.2,-1.5,0.03), 

giving  the  curve  in  Figure  3.6.  With  A  and  C  as  defined  earlier,  the  coefficients  K  are 
given  by 

K  +  A  SC  •  .*012. 
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IV.  ALGORITHM  COMPARISON 


A.  GENERAL 

We  have  mentioned  how  algorithms  are  used  to  fit  the  curve  and  each  algorithm 
has  defects  for  specific  A  values.  Recall  that  the  field  test  program  included  the 
concentrations:  10%,  20%,  40%,  50 %  70%  and  100%  of  full  scale  for  each  element. 
Calibrating  the  tube  consisted  of  two  analyses  each  using  the  20%,  50%  and  100% 
samples;  these  samples  were  used  to  fix  the  absorbance  •  concentration  curve  for  all 
analyses  performed  until  the  next  calibration.  During  the  planned  sequence  of  bums, 
each  of  the  10%,  40%,  and  70%  verification  samples  were  analyzed  twice,  once  just 
after  the  calibration  curve  was  established  and  then  a  second  time  about  30  analyses 
later.  Ideally,  the  established  calibration  curve  should  give  10%,  40%  and  70%  of  full 
scale  respectively  for  these  samples;  the  different  possible  algorithms  can  be  compared 
by  using  these  10%,  40%  and  70%  samples,  in  terms  of  how  close  the  resulting 
calibration  curve  result  is  to  the  known  concentration,  for  the  various  tubes  used. 

Specifically,  each  algorithm  can  be  used  to  establish  a  calibration  curve  using  the 
observed  absorbance  numbers  for  the  20%,  50%  and  100%  samples.  Since  these  curves 
are  not,  in  general,  exactly  the  same  for  all  algorithms,  the  difference  between  the 
returned  concentration  for  a  10%  sample,  e10  say,  can  be  compared  with  the  expected 
result  c,g.  We  will  introduce  the  sum  of  squares  method  to  compare  the  algorithms. 
Thus,  for  a  given  algorithm,  we  will  have  two  e10,  two  e^,  and  two  e70  values  from 
which  we  will  compute  the  sum  of  the  6  squared  differences.  The  algorithm  which  has 
the  smallest  mean  sum  of  squares  will  be  selected  as  the  rank  1,  with  ranks  2,  3,  ....  5 
assigned  the  larger  mean  sum  of  squares.  The  number  of  rank  1  and  rank  2  for  each 
element  throughout  all  tubes  will  be  considered  as  the  comparing  tool  for  that  rank. 
The  algorithm  which  produces  the  most  rank  1  and  rank  2  scores  will  be  designated  the 
best  one. 

B.  SUM  OF  SQUARES  METHOD 

We  will  explain  the  sum  of  squares  method  by  an  example.  For  tube  401,  for 
element  Fe,  the  fitted  curve  is  shown  in  Figure  4. 1  where  the  coefficients  are  computed 
by  algorithm  2. 

The  residuals  E  are  defined  by 
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Figure  4.1  The  Horizontal  Deviation  of  Residual  by  Algorithm  2. 

E  •  C  -  C 

where  C  is  a  known  concentration  and  C'  is  the  estimated  concentration  produced  by 
the  algorithm. 

For  tube  401,  for  Fe,  the  absorbance  data  of  the  verification  sample  A  are, 

A  -  (  28,  105,  129  ) 
and  the  known  concentrations  are, 

C  -  (  20,  50,  100  ) 

Then  we  find  K  by  algorithm  2, 

K  -  (  -0.3211,  0.00602,  0.00159  ) 

We  can  compute  the  estimated  concentration  C'  from  equation  3.4  with  (28,  105,  129) 
in  place  of  A.  That  is 

C  -  (  9.315.44.004,67.040  ). 
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and  the  sum  of  squares  SS  is  defined  as  the  sum  of  the  squares  of  these  residuals, 
i.  eM 

3 

SS  -  £  e-,2  -  0.6852  +  4.0042  +  2.962  -  25.263.  (4.1) 

i»l 


Then  the  mean  sum  of  squares  MSS  is  defined  as, 

MSS  -  (4.2) 

giving 

MSS  -  25.263  /  3  -  8.42. 

Residuals  at  each  point  are  shown  in  Figure  4.1  where  the  residual  is  indicated  by 
the  horizontal  line  between  symbol  '  x '  and  ’o'.  Estimated  concentration  C'  is 
indicated  by  the  symbol  *  x  '  and  concentration  is  indicated  by  the  symbol  'o', 

When  we  are  using  algorithm  4  with  the  same  data,  then 

K  «  (  79.237,-  1.775,0.013  ) 

This  curve  is  shown  in  Figure  4.2.  Horizontal  deviation  is  the  distance  between  symbol 
'  x  '  and  ’o’,  which  is  shown  as  a  horizontal  line.  The  estimated  concentrations  are 

C'  -  (  39.75,  36.57,  67.184  ). 

As  we  see,  the  estimated  concentration  at  absorbance  28  is  larger  than  at  absorbance 
105,  because  the  curve  is  bent  downward  at  the  lower  concentration  with  this 
algorithm.  Therefore,  the  residuals  are  larger.  Also,  algorithm  2  is  considerably  better 
than  algorithm  4  for  this  set  of  data. 

Accordingly, 
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Figure  4.2  The  Horizontal  Deviation  of  Residual  by  Algorithm  4. 

E  -  ( -29.75,  3.43,  2.816  ), 

SS  -  (-29.75)2  +  3.432  +  2.8 162  -  735.257, 
and 

MSS  »  735.257/3  -  245.086. 

For  these  data  with  these  two  algorithms, 

MSS  by  algorithm  2  »  8.42, 

MSS  by  algorithm  4  *  245.086, 

and  algorithm  2  is  better  than  algorithm  4. 

Estimated  concentrations  C'  by  6  algorithms  for  Fe  are  summarized  in  Table  4. 
Known  concentrations  C,  shown  in  the  first  column,  were  used  4  times  10%,  40%  and 
70%  of  full  scale.  After  the  calibration  curve  was  established,  each  of  the  10%,  40%, 
and  70%  verification  sample  was  analyzed  twice.  For  tube  401  this  calibration  was 
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done  twice,  so  that  12  verification  data  are  gathered  in  Table  4.  Notice  the  differences 
between  the  estimated  concentrations  in  column  2  and  3  in  Table  4.  The  estimated 
concentrations  by  PWMA  and  by  algorithm  1  are  very  close.  Since  algorithm  PWMA 
and  algorithm  1  are  exactly  the  same  as  noted  earlier;  the  PWMA  concentration  values 
are  those  printed  out  by  the  microprocessor,  while  the  algorithm  values  are  computed 
from  the  curve  established  from  the  rounded  absorbance  values.  Since  algorithms  2 
through  5  must  be  employed  with  the  rounded  absorbance  values,  rather  than  the 
single  precision  values  available  to  the  microprocessor  the  rank  comparisons  will  be 
restricted  to  only  algorithms  1  through  5.  The  mean  square  residuals  will  be  computed 
for  the  PWMA,  but  they  will  not  be  used  in  comparison  with  the  others;  cases  in  which 
the  PWMA  actually  gives  the  smallest  mean  square  will  be  denoted  by  an  asterisk^). 
We  might  expect  that  algorithm  PWMA  and  algorithm  l  should  perform  very 
similarly.  In  the  last  column,  .00  means  we  can  not  obtain  acceptable  estimated 
concentrations  with  this  algorithm.  If  K22  -  4  K3  (  Kj-  A  )  was  negative,  the  estimated 
concentration  was  forced  to  be  .00  for  all  verification  samples  with  this  calibration. 
The  estimated  concentration  in  the  first  three  rows  of  algorithm  2  and  4  are  the  same 
as  our  previous  example. 

From  the  data  in  Table  4,  we  calculate  the  sum  of  squares  and  the  mean  sum  of 
squares,  which  are  displayed  in  Table  5.  It  says  that  the  sum  of  squares  by  PWMA  is 
100,  and  the  mean  sum  of  squares  is  8.33,  so  the  rank  of  the  mean  sum  of  squares  is 
labeled  with  *  as  it  is  the  smallest.  If  we  don't  count  the  PWMA  in  comparison, 
algorithm  3  is  the  best.  It  has  the  smallest  mean  sum  of  squares,  9.09;  thus  algorithm 
3  suggested  by  Limbeck  is  the  best  for  the  tube  PPM401  for  element  Fe.  The  number 
of  acceptable  data  points  for  each  algorithm  is  given  in  the  row  labeled  N;  if  the  test 
plan  was  followed  exactly,  and  the  algorithm  performed  satisfactorily  for  the 
calibration  samples,  N  would  equal  12  for  each  case.  If  the  algorithm  performed  totally 
unsatisfactorily,  its  N  value  is  0  and  its  mean  sum  of  squares  was  set  equal  to  888  (see 
algorithm  5  in  Table  5). 

The  sum  of  squares,  the  mean  sum  of  squares,  and  the  rank  of  the  algorithm  for 
tube  401  for  9  elements  are  summarized  in  Table  6.  The  equivalent  data  for  the  other 
tubes  are  attached  in  Appendix  D. 

Table  6  show's  the  best  algorithm  for  each  element.  Algorithm  3  is  the  best 
algorithm  for  element  Fe.  Algorithm  1  has  the  smallest  mean  sum  of  squares  for  Ag, 
and  so  on.  As  can  be  seen,  the  best  algorithm  differs  according  to  the  element.  For 
tube  401,  no  single  algorithm  is  the  best  for  all  the  elements. 


TABLE  4 


ESTIMATED  CONCENTRATION  BY  6  ALGORITHMS 


Standard 

Algorithm 

PWMA 

1 

2 

3 

4 

5 

10.00 

8.87 

9.  31 

13.  44 

39.  75 

.00 

E&a 

44.  00 

44.  79 

44.00 

45.  66 

36.  57 

.  00 

EZEH 

cxa 

70.23 

67.04 

67.96 

67.  18 

.  00 

10.00 

13.00 

13.  21 

13.  57 

17.  64 

29.  09 

.  00 

40.  00 

42.00 

42.  32 

41.36 

43.  08 

32.  93 

.  00 

70.00 

70.  00 

70.  23 

67.04 

67.  96 

67.  18 

.  00 

EEB 

8.00 

8.00 

8.  34 

17.04 

53.22 

.00 

Esa 

43.  05 

42.  64 

60.  82 

36.  86 

.  00 

70.  00 

74.  00 

75.  32 

74.  43 

88.  77 

83.  90 

.00 

10.00 

6.00 

6.  07 

6.  39 

13.  58 

66.29 

.00 

40.  00 

36.  00 

35.  69 

35.  19 

52.  89 

26.  16 

.  00 

70.  00 

73.00 

73.  52 

72.  64 

87.41 

81.  58 

.  00 

TABLE  5 

ALGORITHM  COMPARISON  FOR  FE 


Algorithm 

PWMA 

1 

2 

3 

4 

5 

SS(  Fa) 

100. 00 

128.  08 

121.  12 

109. 09 

3119.08 

.  00 

N 

12.00 

Exa 

12.  00 

12.  00 

12.  00 

.  00 

MSS 

8.  33 

10.  67 

10.  09 

9.09 

259.92 

888. 00 

RANK 

* 

3 

2 

1 

4 

5 

In  Chapter  II  we  mentioned  that  Langley  used  6  tubes,  Myrtle  Beach  used  5 
tubes,  Pensacola  used  1 1  tubes,  and  Elmendorf  used  6  tubes.  Due  to  missing  data,  we 
excluded  several  tubes  in  our  analyses.  For  example,  tube  PPM406  doesn't  have 
verification  sample  70  data.  For  similar  reasons  we  discard  1  tube  from  Langley,  1 
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TABLE  6 

SUMMARY  OF  THE  RANKS  FOR  TUBE  401 


Algorithm 

PWMA 

1 

2 

3 

4 

5 

SS( Fe) 

100.  00 

128.  03 

121.  08 

109. 09 

3119. 08 

.00 

n 

12.  00 

12.  00 

12.  00 

12.  00 

12.  00 

.  00 

MSS 

8.  33 

10.  67 

10.  09 

9.  09 

259. 92 

888. 00 

RANK 

* 

3 

2 

1 

4 

5 

SS(Ag) 

2.00 

.  85 

■fl 

12.34 

.  67 

n 

12.  00 

12.  00 

SiKkgiiTtH 

12.  00 

6.  00 

MSS 

.  17 

.  07 

1.  03 

.  11 

RANK 

1 

Hu 

5 

2 

SS(A1) 

268.  00 

317.  88 

300.21 

297.  57 

429. 29 

.  00 

n 

12.00 

12.00 

12.  00 

12.00 

12.  00 

.  00 

MSS 

22.33 

26.  49 

25.  02 

24.  80 

35.  77 

888.  00 

RANK 

* 

3 

2 

1 

4 

5 

SS(  Cr) 

16.00 

17.  29 

17.  16 

17.  10 

39.32 

.  00 

n 

12.00 

12.  00 

12.00 

12.00 

12.00 

.  00 

MSS 

1.  33 

1.44 

1.  43 

1.  42 

3.28 

888.  00 

RANK 

* 

3 

2 

1 

4 

5 

SS(Cu) 

106.  00 

102.05 

100.  06 

99.74 

100. 09 

107. 97 

n 

12.  00 

12.00 

12.00 

12.00 

12.  00 

12.  00 

MSS 

8.  83 

8.50 

8.34 

8.  31 

8.  34 

9.  00 

RANK 

4 

2 

1 

3 

5 

SS(Mg) 

91.00 

170. 50 

156.21 

148. 12 

179. 10 

.  00 

n 

12.  00 

12.00 

12.  00 

12.00 

12.  00 

.  00 

MSS 

7.  58 

14.21 

13.  02 

12.34 

14.  92 

888.  00 

RANK 

* 

3 

2 

1 

4 

5 

SS(Ni) 

38.  00 

34.  87 

26.  02 

31.  70 

192. 41 

.  00 

n 

12.  00 

12.  00 

12.  00 

12.  00 

12.  00 

.  00 

MSS 

3.  17 

2.  91 

2.  17 

2.  64 

16.  03 

888.  00 

RANK 

3 

1 

2 

4 

5 

SS( Si) 

162.00 

149. 09 

141.  29 

166.  29 

171.  89 

103. 04 

n 

12.  00 

12.  00 

12.  00 

12.  00 

12.  00 

6.  00 

MSS 

13.  50 

12.  42 

11.  77 

13.  86 

14.  32 

17.  17 

RANK 

2 

1 

3 

4 

5 

SS(Ti) 

135.00 

142. 77 

364.  12 

141.  42 

134.  41 

141.  47 

n 

12.  00 

12.  00 

12.  00 

12.  00 

12.  00 

12.  00 

MSS 

11.  25 

11.  90 

30.  34 

11.  79 

11.  20 

11.  79 

RANK 

4 

5 

2 

1 

3 

tube  from  Myrtle  Beach,  2  tubes  from  Pensacola,  leaving  24  tubes  in  all.  The  summary 
of  24  tubes  for  each  element  is  shown  in  Table  7  through  Table  15.  We  may  say  that 
algorithms  which  have  the  most  rank  1  or  rank  2  scores  but  the  least  rank  4  or  rank  5 

33 


TABLE  8 

RANKS  OF  24  TUBES  FOR  AG 


scores  are  good  ones,  and  algorithms  which  have  the  least  rank  l  or  rank  2  scores  but 
the  most  rank  4  or  rank  5  scores  are  not  good  ones.  The  entries  in  these  tables  are  the 
counts  of  the  ranks.  For  example,  in  Table  7  algorithm  1  for  Fe  has  7  rank  I  s,  1  rank 
2,  11  rank  3's,  5  rank  4's,  and  0  rank  5  out  of  24.  For  Fe,  the  PWMA  internal 
algorithm,  with  rounded  concentrations,  had  the  smallest  mean  sum  of  squares  10  time 
out  of  24.  In  Table  7,  we  conclude  that  algorithm  1  is  better  than  algorithms  2,  3,  4, 
and  5.  Algorithm  2  has  7  rank  l's  but  also  5  rank  5's.  Therefore,  algorithm  2  is  not  as 
good  as  algorithm  1.  If  algorithm  2  did  not  have  those  undesirable  cases  with 
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TABLE  9 

RANKS  OF  24  TUBES  FOR  AL 


denominator  zero,  then  it  would  be  the  best.  Also  algorithm  4  has  11  rank  5's  because 
it  has  a  lot  of  deviations  for  two  extreme  cases  as  we  saw  in  the  previous  Chapter. 
Algorithm  5  is  almost  the  same  as  algorithm  4.  Algorithm  3,  with  5  rank  I  s  and  12 
rank  2's,  is  acceptable  but  not  as  good  as  algorithm  1.  We  conclude  that  algorithm  1 
is  the  best  for  element  Fe. 

For  element  Ag,  even  though  algorithm  5  has  9  rank  I  s,  it  also  has  6  rank  5's. 
Therefore,  we  can  not  select  algorithm  5  as  the  best  one.  Algorithm  2,  or  1  are 
acceptable  for  element  Ag. 
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For  Al,  algorithm  2  dominates  the  other  algorithms,  since  it  has  10  rank  I  s. 
Algorithm  3  is  the  second  algorithm.  Algorithm  3  could  be  useable  even  though  it  has 
fewer  rank  1  scores. 

For  Cr,  algorithm  4  is  not  acceptable,  for  it  has  22  tubes  rank  4's  or  5's. 
Algorithm  1  is  the  best  for  element  Cr. 

For  Cu,  Mg,  and  N'i,  algorithm  1  dominates  the  others.  For  Si,  algorithm  2  or  3 
are  the  best.  For  Ti,  algorithm  2  is  good.  Algorithm  1  and  4  are  acceptable  too. 


TABLE  13 

RANKS  OF  24  TUBES  FOR  M 


Algo. 

rank 

* 

1 

2 

3 

4 

5 

PWMA 

7 

1 

13 

2 

8 

1 

0 

2 

4 

2 

2 

5 

11 

3 

4 

18 

2 

0 

0 

4 

1 

1 

3 

16 

3 

5 

2 

1 

9 

2 

10 

TABLE  14 

RANKS  OF  24  TUBES  FOR  SI 


Algo. 

rank 

* 

1 

2 

3 

4 

5 

PWMA 

11 

1 

3 

7 

8 

3 

2 

2 

8 

7 

2 

3 

4 

3 

7 

2 

3 

11 

1 

4 

1 

7 

9 

4 

3 

5 

4 

1 

2 

3 

13 

Table  16  summarizes  the  results  of  the  previous  tables  by  listing  the  best  and 
alternate  algorithm  for  each  element.  We  can  say  that  algorithm  1  provides  the  best  fit 
for  elements  Fe,  Cr,  Cu,  Mg,  and  Ni,  and  algorithm  2  provides  the  best  fit  for  elements 
Ag,  Al,  Si,  and  Ti. 
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TABLE  15 

RANKS  OF  24  TUBES  FOR  TI 


Algo. 

rank 

* 

1 

2 

3 

4 

5 

PWMA 

11 

1 

5 

9 

3 

5 

2 

2 

7 

3 

2 

8 

4 

3 

4 

1 

6 

7 

6 

4 

5 

5 

11 

2 

1 

5 

3 

5 

3 

2 

11 

TABLE  16 

THE  BEST  ALGORITHM  FOR  EACH  ELEMENT 


Element 

Fe 

Ag 

A1 

Cr 

Cu 

Mg 

Ni 

Si 

Ti 

First 

1 

2 

2 

1 

1 

1 

1 

2 

2 

Second 

3 

1 

3 

3 

2 

3 

3 

3 

n 

V.  CONCLUSIONS  AND  RECOMMENDATIONS 


There  is  no  single  best  algorithm  for  all  nine  elements,  since  the  degree  of 
curvature  is  different  according  to  the  element.  The  PWMA  didn't  consider  different 
algorithms  for  different  elements.  The  computation  in  the  PWMA  used  only  by  one 
algorithm  for  all  nine  elements.  If  we  employ  different  algorithms  for  different 
elements  in  an  analyzer,  it  will  be  better  because  each  element  has  different 
characterisics.  We  recommend  that  one  may  use  algorithm  1  for  elements  Fe,  Cr,  Cu, 
Mg,  and  Ni  and  algorithm  2  for  elements  Ag,  Al,  Si,  and  Ti  to  obtain  better  estimation 
in  computing  concentration.  The  microprocessor  can  easily  check  for  the  existence  of 
an  asymptote.  If  one  occurs,  use  algorithm  1. 


APPENDIX  A 

PWMA  HELD  TEST  BURN  SEQUENCE 


Iky  1 

1.  Tub*  Condition 

2.  Auto  Zoro  -  Op* rata  in  A Z  aode.  If  th*  AZ  ftiladi  rap oat  tha  A Z  and 

nudaar  tha  run  as  ZA.  If  tha  AZ  failad  twice  in  a  row, 
replace  tha  tuba. 

S.  CAL  la  -  Operate  in  CAL  X  aode. 

4.  CAL  Zb  -  Operate  in  CAL  I  node. 

5.  CAL  lie  -  Operate  in  CAL  IX  aode. 

4.  CAL  lib  -  Operate  in  CAL  IX  aode. 

7.  CAL  Ilia  -  Operate  in  CAL  III  node. 

8.  CAL  ZlXb  -  Operate  in  CAL  III  aode. If  tha  calibration  failad, 

repeat  tha  calibration  and  nuaber  tha  calibration  run*  as 
3A,  8A,  accordingly. 

9.  Blank  -  Run  a  bum  with  no  saaplai  operate  in  Analysis  node.  Fro* 

now  on,  all  tha  aaaples  will  be  run  in  Analysis  node  except 
Res lop*. 

10.  Verification  Saapla  8  1:  10X 

11.  Verification  Saapla  8  2:  40Z 

12.  Verification  Saapla  8  3s  70Z 
13  Blank 

14- 18. Correlation  saapla  8  1  -  Repeat  tha  saaa  saapla  five  tiaas,  this 

procedure  applies  for  all  tha  oor relation  saapla. 

19-23. Correlation  saapla  8  2 
24-28. Correlation  saapla  8  3 
29.  Blank 

30-37. Randoa  saapla*  -  These  aay  be  real  angina  oil  saapla.  If  there  are 

not  enough  real  angina  oil  saapla* ,  repeat  one  saapla 
as  aany  tiaas  as  necessary.  (  8  saapla*  total  ) 

38.  10X  saapla 

39.  40Z  saapla 

40  70Z  saapla 

Day  2 

41*  Blank 

42.  Raslopa  A  -  Operate  in  RSLP  aode.  If  tha  res lope  A  failad,  disregard 

and  continue  tha  sequence. 

43.  Raslopa  8  -  Operate  in  RSLP  aode.  If  th*  raslopa  8  failad,  disregard 

and  oontinue  tha  sequence. 

44.  10Z  saapl* 

45.  40Z  saapla 

44 .  70Z  saapla 

47- 51. Correlation  saapla  8  4 
52-54. Correlation  saapl*  8  5 
57.  Blank 

58-43. Randoa  saapla*  (  4  saapla*  total  ) 

44.  CAL  I  chock  saapl*  -  Operate  in  Analysis  aode. 

45.  CAL  II  check  saapla  -  Operate  in  Analysis  aode. 

44.  CAL  III  check  saapla  -  Operate  in  Analysis  nods. 

47.  Blank 

48- 77. Randoa  saapla*  (  10  saapla*  total  ) 

78.  10Z  saapla 

79.  40Z  saapl* 

80.  70Z  saapla 

Oay  3 

81.  Blank 

82.  CAL  la. 

83.  CAL  lb. 
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B4. 

CAL  ZZa. 

85. 

CAL  ZZb. 

B4. 

CAL  ZZZa. 

B7. 

CAL  ZZZb. 

80. 

Blank 

09. 

102  aaaple 

90. 

402  aaepla 

91. 

702  aaaple 

*2-**.  Correlation  Map  la  0  4 
97-101 .  Correlation  i«pU  0  7 
102.  Black 

103-117. Randoe  taplM  <  U  aaaplea  total  ) 
110.  102  Maple 

119.  402  saapla 

120  .  702  Maple 

Day  4 

121.  Black 

122.  Realope  A 

123.  Realope  B 

124.  102  Maple 

125.  402  aaepla 

124.  702  Maple 

127-131. Correlation  aaaple  0  B 
132-134. Correlation  aaaple  8  9 
137.  Black 

138-1B7.  Randoe  aaaple*  (  20  saapla*  total  J 
ISO.  CAL  Z  check  aaaple 

159.  CAL  ZZ  check  aaaple 

160.  CAL  1ZZ  check  aaaple 
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APPENDIX  B 
DATA  PPM401 


There  are  several  data,  we  list  just  one  of  them  for  references. 

PPM401 


ms  4oi 

CAL  IA 

9 

S 

00 

OS 

10 

58 

25 

44 

40 

10 

10 

cal  n 

9 

4 

04 

02 

10 

54 

22 

44 

58 

8 

4 

CAL  XXA 

9 

0 

110 

48 

14 

85 

47 

40 

42 

28 

50 

CAL  XXO 

9 

4 

111 

100 

21 

84 

45 

41 

44 

18 

24 

CAL  XXXA 

9 

7 

145 

124 

24 

104 

70 

72 

114 

55 

42 

cal  xxxa 

9 

0 

140 

151 

24 

104 

71 

74 

120 

55 

44 

VCRIF  10 

9 

10 

9 

1 

2 

1 

4 

5 

5 

2 

5 

VCRXF  10 

AOSOROANCCS 

28 

28 

0 

17 

12 

52 

14 

4 

11 

VCRIF  40 

9 

11 

44 

4 

10 

4 

18 

4 

12 

4 

7 

VCRXF  40 

AOSOROANCCS 

100 

40 

17 

70 

42 

04 

74 

20 

20 

VCRXF  70 

P 

1C 

70 

7 

.  4 

4 

28 

17 

20 

12 

10 

VCRXF  70 

AOSOROANCCS 

124 

114 

14 

88 

04 

44 

105 

24 

24 

VCRXF  10 

P 

so 

15 

1 

8 

1 

8 

2 

4 

4 

8 

VCRXF  10 

AOSOROANCCS 

40 

28 

10 

20 

22 

50 

25 

21 

21 

VCRXF  40 

P 

39 

42 

4 

7 

0 

14 

8 

12 

15 

7 

VCRXF  40 

AOSOROANCCS 

101 

84 

14 

77 

54 

54 

74 

50 

18 

VCRXF  70 

P 

40 

70 

7 

11 

7 

24 

14 

21 

12 

12 

VCRXF  70 

AOSOROANCCS 

124 

110 

18 

44 

04 

48 

107 

24 

54 

P524  401 

RCSLOK  1 

P 

4C 

117 

101 

12 

72 

01 

55 

80 

57 

14 

RCSLOPC  2 

P 

45 

114 

44 

15 

44 

01 

54 

80 

54 

17 

VC RIF  10 

P 

44 

12 

1 

1 

1 

4 

2 

5 

4 

5 

VCRXF  10 

AOSOROANCCS 

54 

24 

5 

17 

15 

28 

14 

10 

7 

VCRXF  40 

SX>LXNXT  P 

40 

44 

4 

8 

4 

20 

8 

11 

25 

8 

VCRXF  40 

AOSOROANCCS 

107 

41 

10 

74 

40 

55 

70 

42 

21 

VCRXF  70 

SX>LIHXT  P 

44 

74 

4 

10 

7 

51 

14 

17 

54 

10 

VCRXF  70 

AOSOROANCCS 

151 

124 

14 

40 

41 

45 

44 

47 

50 

CAL  X 

P 

44 

25 

2 

5 

2 

8 

4 

0 

8 

4 

CAL  X 

AOSOROANCCS 

45 

02 

7 

58 

25 

42 

54 

17 

10 

CAL  IX 

P 

40 

00 

0 

4 

4 

20 

14 

14 

14 

4 

CAL  XX 

AOSOROANCCS 

115 

48 

14 

88 

44 

41 

84 

28 

25 

CAL  XXX 

CR  »NX,SX»P 

44 

105 

10 

28 

11 

50 

25 

51 

50 

25 

CAL  XXX 

AOSOROANCCS 

144 

150 

28 

110 

45 

74 

124 

50 

50 

VCRXF  10 

SI>LXHIT?P 

78 

4 

1 

2 

1 

4 

2 

2 

0 

s 

VCRXF  10 

AOSOROANCCS 

28 

24 

0 

21 

12 

50 

15 

117 

4 

VCRIF  40 

P 

79 

41 

4 

11 

0 

15 

11 

12 

17 

8 

VCRXF  40 

AOSOROANCCS 

44 

80 

18 

74 

57 

54 

74 

51 

25 

VCRIF  70 

SI>LINXT  P 

80 

75 

7 

22 

7 

20 

25 

21 

25 

15 

VCRIF  70 

AOSOROANCCS 

151 

114 

24 

40 

55 

72 

107 

3.9 

54 

APPENDIX  C 
PROGRAMMING  LIST 


7 COEP1 

fl  8/12/86  PROGRAM  TO  FIT  BILL « S  2  POINT  RATIONAL  EQUATION.  CHECK  OUT 
A  VERIFICATION  READINGS  AND  THE  VALUES  AFTER  THE  RESLOPE  ks  NELL . 

H+  10. 10. 25  0.10. 40. 250. 30. 20. 2 
#4-20 +5  xx  9x1+1 
T1+  8  3  *  6  1  pO 
R1+-  3  0  *  6  1  pO 

914*6  0  *  6  1  pO 
FOR  TUBE  401 « 

'MODEL  A/ IK1+K2A)' 

LOOPIA+W3  2p(*,PPM40lCl  +  x6j£CI]+t5]  )  +  2 
□♦JftCULl  2]*  20  50  xffCijJiiCl  2j  *  .  *  0  1),[0.1]  (AC2  33*  50  100 
■flCJJ 5  All  3]«.*  0  1 

a  ORIGINAL  COEFFICIENTS  FOR  THE  TWO  CURVES. 

D*-#  .PPM401  L7+xl2  j£[I]+x5 j 
©♦PC2xi6]*(,(PC2xi6j«».*  0  l)  +  .x$JOCl  3  6  7  9  12] 

U1+  8  3*61  p0_ 

Q+6  0  *  6  1  oDZ  1+2xi6] 

01+QlfCl]§ 

04.(6  0*6  lpflC  l+2xX6]),01 

ri4-ri,  ci]  ui 

P4-  3  0  *  6  1  p  10  40  70  Xff  C  (PCI]  *5  )-4] 

P14-P1,  rijfl 

'FOR  ELEMENT: » .  (ELlIil  Y.\.  PRIOR  TO  RE  SLOPE ' 

U+U+C3  0  *  6  1  p  10  40  70  x&t(feci]+5 )-4] ) , 0 
C* + /  (  * .  PPM4  0  i  C  2  0  + 1 2  s  E  [  I ]  + 1 5 1 )  ♦  2 
C4.50xflCn  +  (2pCJ*(X+.xC*  0  1) 

A  GIVES  MULTIPLICATIVE  CORRECTION  FACTORS  FOR  THE  TWO  LINES 
P4-#  PPAI401 C (22+1 6  J . 34+t 6 iPCi’] +i 5] 

^£££6  CjD(2xt6]  °  0  1  )+.x4jK)[l  3  6  7  9  12] 

•  AFTER  THE  RESL.OPE ' 

CM6  0*61  oDZ  l+2xi6]  ).  8  3  *  6  1  pO 
D+* .PPM401 L (28+i 6 )  jffill  +  x  5] 


'1+2x16] ),U 1 


U+U+{3  0  *  101  p  10  40  70  20  50  50  100  10  40  70  xfl£(, 
+.((I4-I+1)S9)/£00P 
1 4-1 

LOOTzA 4-+/  3  2  p(«,PPM401j:41+x6;FCI]+\5]  )  +  2 
K+CCACl  2]+2BWxffCI]  JlACl  2]  •.*  0  l),t6.1]  Ul 
fflffLJU)  All  3]  •  .  *  0  1 
04-*  .PPM401 C47  +  112  jPCJ]  +  x  5] 

0+PC2xi6]  +  (,(PC2xi6]o.*  0  l)  +  .xW)Cl  3  6  7  9  12] 
U 14-  8  3*61  P0_ 

Q4-6O  *6  I9PC  1+2xi6] 

Ql+Ql , Cl]  Q 

u4-(6  0*61  p0C  l+2xx6]  )  tUl 
T1+T1,111U1 

P4-3  0  »  6  1  p  10  40  70  xffC(PCJ]  +  5)-4] 

Rl+Rl.  M.]  R 

•  FOR  ELEMENT: ' .  {EL  II:  ]  ) ,  >  .PRIOR  TO  RESLOPE ' 
□+04.f3  0  *  6  1  p  10  40  70  xflt(feci]+5)-4] ),0 
4.((l4-I+l)S9)/L00r 

7  COEFR 


(AC2  3] *  50  100 


8/12/86  PROGRAM  TO  FIT  BILL '  5  2  POINT  RATIONAL  EQUATION .  CHECK  OUT 
VERIFICATION  READINGS  AND  THE  VALUES  AFTER  THE  RESLOPE  ks  WELL . 


H-*-  1  0.1  0.25  0.1  0.4  0.25  0,3  0.2  0.2 

E4-20+5xt9xI«-l 

T24-  8  3  *  6  1  pO 

'FOR  TUBE  401 ' 


IFO^  THE  TWO  CURVES . 

3|5xCC2x  X6]x0[2xx6]  )+((EC2 j]x0[2xx6]  )-l) 


oT  7  It •  £  ljoi  "1+2xx6]  ),£/2 

O+/U.PPH40ip0  +  x2:E[I]+x5y)  +  2 

05  0  xff  f  n  ♦7m j  ]  xc+it  3 : 1  xcxc  j  ♦  ?  (X[  2 1 T xC 0  - 1 ) 

A  CIVES  MULTIPLICATIVE  c6rRECTI0N  FACTORS  FOR  Tl 


? SLOPE' 


]  A  GIVES  MULTIPLICATIVE  CORRECT 
I  O* .PPM401  C(22  +  x6  5.34+x6|E£i’ 
_ J  OCL60I]  x  (aLI  i  ]  xi?cixi6  J  +aL3  t 

;  2  3  •  AFTER  TEE  RESLOPE  ' 

24  O(7  1*6  1p0[  1+2x16]).  8  3  « 

OA.PPM40U(28+x6):E[.n+x5] 
OCC4pl]x(X[ij]xnfi  44  6]+tf[2 


fjiici 


'ORS  FOR  THE  TWO  LINES 
6]xDC2xi6])+((irC2j]x£)C2 


1+2x16]).  8  3*61  p 0 
+i6)iECI]+\5] 

]xDC2  4  4  6]+tfL3i]x0[ 


,25]  O«,PPM40lt(28+x6)iECI]+\5] 

:26]  OCC4pl]xUCljJxPvt2  4  4  6]+tfC3l3x0[2  4  4  6]x0[2  446]  )♦((*[ 
?!J  .PC?  4  4  6]  )-l) 


♦  l  ( t/022  000  f=  1 )  NRR 

♦Tail +1^9  f/loop*^  ^0C503100  1  0*4  0  7  0  xA^fJ  +  J  f-4  ]V,U 


;32]  LOO 
33]  K*( 


(• #PPM401 C4l+l6  jPCJ] +l 5] )+2 
,A*2)tX 2i)  50  100  xfftJl  )o.*  1  0 


)+, t^2i)  50  100  xfftJ]  )•.*  1  0  1 
EFFICIENTS  FOR  TEE  TWO  CURVES . 


a  ORIGINAL  COEFFICIENTS  FOR  THE  Tl 
O*  .PPM40 1  £47  + 1 1 2 1 El f]  + 1 5  J 
0(atlj]xfiC2xi6])+XC3|]x5C2xx 
(72+  8  3  *  6  1  p0  „ 


6]xD[2xx6]  )+((JfC2j]xDC2xx6]  )-l) 


O  (71*61  pD  C 
T2*T2jlZ11  U2 


l+2xx6] ) ,02 


WELL . 


'FOR  ELEMENT : '  ,  (ELCIx]  ) PRIOR  TO  RESLOPE ' 

7i  «*  <«*•»>-«).* 

;44,  -*-0 

,45  E/?E:JJ+4p0  ,  , 

.46.  0(  3  0  +0),[1](<<7  1  *  4  1  pPCl  3  3  5]  ),  8  3  *4  1  00). [1]  3  0  +0) 
0+0(4  1  *  10  1  p  10  40  70  20  50  50  100  10  40  70  xff[  (ELlJ  +  5  )-4]  )  ,U 
♦( (I+J+l )*9 )/L00P 

V  C0EF3 

A  9 / 0 6 /  8 6  PROGRAM  TO  FIT  BILL'S  2  POINT  RATIONAL  EQUATION .  CHECK 
A  VERIFICATION  READINGS  AND  THE  VALUES  AFTER  THE  RBSLOPE  AS  WELL 
A  MODEL  C  =  A+(Cl+C2xil+C3xAxA  ) 

'MODEL  C-  A*  lci+C2*A+C3xAxA  )  * 

E+  10. 10. 25  0.10. 40. 250. 30. 20. 2 
T  3+  8  3*6  1  p0 
P+20+5xx9xJ+l 
'FOR  TUBE  401 » 

LQOPSA++/  3  2  p  (*  .PPM201  [I  +  16 }E[J] +x  5]  )+2 
O+E+U+O  20  50  100  xff  rmifto  .  *  0  1  2 
A  ORIGINAL  COEFFICIENTS  FOR  THE  TWO  CURVES . 

Oa .PPM401 l7+x 12 ;PCr]+x  5] 

0+PL2xx6]  +  (,  (PC2xx6]  0  .★  0  1  2  )+.x$E) 

U 3+  8  3  *  6  1  pO  _ 

0(60*6  lpLC  l+2xx6]  )  ,U3 
T3+T3. Cl]  U3 

'  FOR  ELEMENT: '  .  (EL  [Is  ]  )  .  •  .PRIOR  TO  RESLOPE ' 

□+0(3  0  *  6  1  p  10  40  76  x&V(E[Il  +  5)-4]  ),(/ 

0+/ U. PPM401 C20  +  X2  }E[I] +t  5]  )  +  2 
C+50+(lpC)+(E+.xC*  0 1 2 J 

A  GIVES  MULTIPLICATIVE  CORRECTION  FACTORS  FOR  THE  TWO  LINES 
Oa.PPM401C(22  +  x6),34+x6  jECl]  +  x5]  v 
OCC6pl]xZJC2xx6]  +  (;(LC2xx6]o.*  0  1  2)+.x$K) 


'AFTER  TEE  REi 


U*(  6  0*  6  1  oDt  1+2x163  )4  8  3  16  1 
D* 1 .  PPm  0 1  [  (  2  8  + 1 6  )j  £  [12 + x  5  ] 

WlhWMtioV0  70  2a  50  50  1 


,  (DC2x  1  2  2  32  8.*  0  12  +.xW 
>  1  POCI  3  3  5]  ,  8  3  9  4  1  P0),[12  3  0  +0) 
20  50  50  100  10  40  70  xffHPtrJ  +  S  )-42  )  ,0 


1*1 

LOOT i  A**/  3  2  p  (*  ,PPAf201 
U*K*Ca*C*  20  50  100  xff  [ I 
n  ORIGINAL  COEFFICIENTS 
D** .PPM 401 [47+ 1 12 ; ECJ] 
0+-0L2xi62  +  (,  (0l2xi628. 


+i6{0[I2+i52  )+2 
A®  .*  0  1  2 
TEE  TWO  CURVES . 


0«-0C2xi6]  +  (,O)C2xi6Jo.*  0  12)+.xW) 

03+  8  3  1  6  1  p0  _ 

^  6  0  9  6  lp0[  1+2x162), 03 
T3*T3.cii  03 

'  FOR  ELEMENT : '  .  (ELlIil  ).' .PRIOR  TO  RESLOPE  ' 
U*U* (3  0  *6  1  p  10  40  70  xffI<fe[J]  +  5)-42l,0 
*(  (1*1 +1H9)Y LOOT 


V  COEFU 

n  9/09/86  PROGRAM  TO  FIT  BILL '  S  2  POINT  RATIONAL  EQUATION 
p  VERIFICATION  READINGS  AND  THE  VALUES  AFTER  THE  KESLOPE , 
p  MODEL  C  -  01+02x4+03x4x4 
'MODEL  0  =  01+02x4+03x4x4  * 

H*  1  0.1  0.25  0.1  0.4  0.25  0.3  0.2  0.2 
^♦20+5xigxI^l 


HECK  OUT 
WELL. 


T4+-  8  3  1  6  1  pO 
'FOR  TUBE  401 »  > 

LOOP: A**/ 3  2 o(l.PPM50 
Q*K*XC*  20  50  100  X0JTJ1 
p  ORIGINAL  COEFFICIENT 
D**  .PPM40 1  Jf 7  + 1 1 2 t E tl 

g*,  v0[2xi62  0  .  *  0  1  2)  + 
4+1  3  96  1  00 
T4+r4,ci]04  _ 


L+t6*0[/2+t52 )+2 

4  8  .  *  0  1  2 

OR  THE  TWO  CURVES . 


r4«+r4,ci] 


Lp0C  1+2xi62),04 
.wnJNTi '  .  (ELLlil )  . '  .  P0I< 
3  0  1  6  1  p  lo  40  70  xflWfeC1 14 
l,PPM40i[20+i2;E[I2+i52) 
x0[I2 )♦(£+. x0*  0  1  2J 
5  MULTIPLICATIVE  CORRECT K 
PM4Ol[(22+t6),34+i6i0[I2 
2xi6J  • .  *  0  1  2  )  +  .  xif ) 


mimr 


U*U* ( 3  0  1  6  1  p  10  40  70  xJJT  (2JCI2  +  5  )-4]  )  0 
0+-+/Cl,PPM4Oi[2O  +  i2;E[I2+i5])  +  2 
0+-(5Oxffri2J+(fc+.x0*  0  1  2 j 

o  GIVES  MULTIPLICATIVE  CORRECTION  FACTORS  FOR  THE  TWO  LINES 
0+-1  ,PPM401  [  (22+t6  )  ,34+i6  {0LT2  +t  52 
0*Cx(,(D C2*i 6]  •.*  0  1  2  )  +  . xjf ) 

•  AFTER  THE  RESLPPE ' 

U*( 6  0  16  1  qDI  1+2xi62  ).  8  3  1  6  1  pO 
0+-1  ,PPM401[ (28+i6j  {ElIJ+i  52 
0*0*  ( , (0[2x  1  22  3J8.*  0 1 2  )  +  .  x)Q 

0+-(  3  0  4-0), [12  (((6  0  14  1  p0[l  33  52  ),  83141  p0).[l2  3  0  +0) 
□+-0+-(3  0  1  10  1  p  10  40  70  20  50  50  100  10  40  70  xfl[  (0[I2+5  )-42  )  ,0 


0<-0x(,(0[  2x  1  22  3  j  8  .  *  0  1  2)+.xX) 

U*\  3  0  4-0), [12  (((60141  p0[l  33  52  ),  83141  P0).[12  3  0  +0) 

2(  (lil+Wg0)1  loop**  70  20  50  50  100  10  40  70  xflC(^**  +  5)-4]  )»u 
!♦! 

LOOT i A**/  3  2  p (# .PPM501 [41+t6 J0CI2 +t 52 )  +  2 
!>)C«-(0>*-  20  50  100  x0CI]  Ji48 .  *  0  12 
p  ORIGINAL  COEFFICIENTS  FOR  THE  TWO  CURVES. 

0+-l»PPM4OlL47  +  \12  }0 [12  +  l 52 
0*,  (0 [2xi 62°.*  0  1  2 J  +  .x/f 
04+-  8  3  1  6  1  P0 
T4+T4  ,  [  1 2  04  _ 

0+(6Ol6  1pD[  l+2xa32), 04 
'  FOR  ELEMENT  l '  .  (EL  ill 2  )  . '  .PRIOR  TO  RESLOPE' 

U*U* (3  0  16  1  p  10  40  70  xNt(fe[I2*5)-42l,0 
♦  (  (I-+I+1  ;£9  )/L00T 

7  COEF5 

P  9/09/86  PROGRAM  TO  FIT  BILL '  S  2  POINT  RATIONAL  EQUATION .  CHECK  OUT 
P  VERIFICATION  READINGS  AND  THE  VALUES  AFTER  THE  RESLOPE  AS  WELL . 

P  MODEL  A  -  01+02x0+03x0x0 
'MODEL  A  -  01+02x0+03x0x01 


►  10. 1  0. 25  0.10. 40. 25  0.30. 20. 2 
►20+5*i9*X<*l 


11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
.23. 

24 

25 

26 

.27, 

28 
if  Sc 

.29, 

30 

31 

32 

llti 

.33, 

34 

35 

l”ZI 

36 

■  :zi 

.37, 

38 

39 
>0 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

Iso; 

60 

61 

62 
iii  i 

.63, 

64 

65 

•  %  X  > 

66 

i  X  X  i 

.67, 

68 

69 

iSzi 

.70, 

71 

72 

73 


PPMHO 
00  “ 
cm 

\E  _ 

.<100* _ 

23*2;-4xJf 


+  i6;P[X3+i53)+2 
•.*0  12 

i  THE  TWO  CURVES. 


ami-. 


0[2*i63 ))*0.5)+2*P[33 


■OR^SLi 


SLOPE ' 


?5+ 

T5tr5._ 

U+  {  6  0*6 
'FOR  BLEW 

□♦0^(3  0  4  6  _  r  . 

0+/(f.PPJ?40i[20+i2.  __ 

C+C*lS+.*(50*SCJ]J*  0  12. 
n  GIVES  MULTIPLICATIVE  CORRECTION  FACTORS  FOR  TEE  TWO  LINES 
0J» jPPU401  c  <22+l6  )  ,  34+16  ;P[X3  +l 5 1 

-  )-4xP[33*(P[lD -L[2xXX3  ) 

0*0 . 5 )+2*P[33 

DOQ 

8  3  4  6  1  p0 

+1  53 


i+2*i6jj^ 


8+16 ) jP! 


(/♦< 6  0  4  6  l  ofll 
©♦•jPPjwoitTi 

KK+  2  4  4  6 

__  - *2)-4xP[33*<K[13-0UMCXX33> 


001 : DD+ 
♦  (00< 0 
Oil!}*! 


liUhn+ri.rs 

. . .  . 

MiSiW0  70  20 

J>1 
LOOT : 

□♦*♦ 


0*0 . 5  )+2*P[33 
n01 
23*2 
6  0  4  4 


0  50  50  100  10  40  70  x0[XP[X3+5  )-43  )  ,0 


P[X3+l53)  +  2 
2 


?S  FOR  THE  TWO  CURVES . 

-4llc3lxU?l3-D?lx?6??? 


PI 

*0.5)+2*K[33 


U+{ 6  0  4  6  1  pDl 
jR  ELEMENT  i' 


'Fdi  f 

MtCto^e 


1+2*133  ),  8  3  4 


,  _ _ JRIOR  TO  RES  Li 

„  .  _  1  p  10  40  70  x9I  (P[I3  +  5 )-43 1  ,U 

(U+I+l)S9)/LOOT 

0TBERU05+  8  3  *6  1  pO 
0+<6O46  1p0[  l+2*\  63  ) ,  8  3  4  6  1  pO 
r5«-r5.  [1305 


SLOPE ' 


5^5,  [13  05 

POP  ELEHENTi «  .  (PL  [1: 3  h,'.  PRIOR  TO  RE  SLOPE 1 
k0<-(3  0  4  6  1  p  10  40  70  *»t(P[I3+5)-43  ),0 


;s9S5l00x 

.113 -*8  8  88 


□♦0+(3  0  4 

(  (JVJ+lJ: 

PPP0P1J0CJ 
♦LI 

ERRORQ i 0 CXX3 ♦8888 
♦L0 

«PPP0P1 :0[XX3^8888 
fi^Ll 

fl ERRORQ : 0 [XX3  +8  8  8  8 
fl+LQ 

OTHER: J+I 

0^v6  0  4  6  1  p0[  1+2*163  ),  8  3  4  6  1 
05+  8  3  4  6  1  pO 


P0 


46 


APPENDIX  D 
SUMMARY  OF  RANK 


OATAAL  PPM201 


FE 

175.00 

190.49 

A5.59 

119.33 

972.52 

00 

SIZE 

A. 00 

A. 00 

A. 00 

A. 00 

A.  00 

.00 

RATIO 

29.17 

31.75 

10.93 

19.89 

142.09 

888.00 

RANK 

3 

1 

2 

4 

5 

AC 

20.00 

22.11 

33.74 

22.85 

37.23 

.00 

SIZE 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

.00 

RATIO 

5.35 

3.A9 

5.A2 

3.81 

A. 20 

888.00 

RANK 

II 

1 

3 

2 

4 

5 

AL 

3A.00 

♦9.29 

47.87 

47. 7A 

47.82 

48.75 

SIZE 

A.  00 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

RATIO 

A. 00 

8.22 

7.98 

7.9A 

7.97 

8.12 

RAMC 

» 

5 

3 

1 

2 

4 

CR 

10.00 

10.04 

9. 44 

9.58 

21.59 

.00 

SIZE 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

.00 

RATIO 

1.A7 

1.A7 

1.58 

1.A0 

3. AO 

888.00 

RANK 

3 

1 

2 

4 

5 

CU 

28.00 

29.44 

33.00 

31. 8A 

42.77 

33.90 

SIZE 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

RATIO 

4.A7 

4.91 

5.50 

5.31 

7.13 

5.A5 

RANK 

II 

1 

3 

2 

5 

4 

MB 

32.00 

18.54 

20. 94 

21.13 

117.57 

21.28 

SIZE 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

RATIO 

5.33 

3.09 

3.49 

3.52 

19.59 

3.55 

RANK 

1 

2 

3 

5 

4 

NI 

19.00 

25.  A2 

370. 9A 

14. 47 

89.70 

.00 

SIZE 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

.00 

RATIO 

3.17 

4.27 

A1.83 

2.44 

14.95 

888.00 

RAMC 

2 

4 

1 

3 

5 

SI 

283.00 

.00 

383.29 

193.42 

.00 

.00 

SIZE 

A. 00 

.00 

A. 00 

A. 00 

.30 

.00 

RATIO 

♦7.17 

888.00 

A3. 88 

32.24 

884.00 

888.00 

RANK 

3 

2 

1 

4 

5 

TI 

22.00 

21.31 

1A7.05 

18.08 

21.07 

19.45 

SIZE 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

RATIO 

3.A7 

3.55 

27.84 

3.01 

3.51 

3.24 

RAMC 

4 

5 

1 

3 

2 

OATAAL 
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888.00 

888.00 

RAW 

2 

1 

3 

4 

5 

CR 

2.00 

1.87 

2.05 

1.91 

3.59 

2.19 

SIZE 

6.00 

6.00 

6.00 

6.00 

6.00 

6.00 

RATIO 

.33 

.31 

.34 

.32 

.60 

.36 

RAW 

1 

3 

2 

5 

4 

53 


cu 

IAS. 00 

179. OA 

17A.91 

179.00 

212.23 

17A.33 

sza 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

A.  00 

RATIO 

27. SO 

29.90 

29.49 

29.05 

35.37 

29.39 

RANK 

* 

4 

2 

3 

5 

1 

MO 

110.00 

110. OA 

119. OA 

119.04 

4A7.95 

.00 

sza 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

.00 

RATIO 

19.67 

19.  AO 

19.90 

19.97 

77.99 

000.00 

RAM 

II 

1 

3 

2 

4 

5 

NZ 

3S.OO 

32.50 

27.04 

21.01 

01.19 

.00 

sia 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

.00 

RATIO 

5.03 

5.42 

4.A4 

3.A4 

13.53 

000.00 

RANK 

3 

2 

1 

4 

5 

SZ 

19.00 

14.10 

17.20 

1A.45 

19.92 

1A.95 

sza 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

RATIO 

3.17 

2.3A 

2.00 

2.74 

3.32 

2.03 

RANK 

1 

4 

2 

5 

3 

TI 

30.00 

51.00 

1A07.75 

A1.A5 

51.55 

.00 

SIZE 

A. 00 

A. 00 

A. 00 

A. 00 

A. 00 

.00 

RATIO 

A. 33 

0.51 

201.29 

10.20 

0.59 

000.00 

RAM 

K 

1 

4 

3 

2 

5 

DATAAL 

PPM401 

FE 

100.00 

120.03 

121.00 

109.09 

3119.00 

.00 

sia 

12.00 

12.00 

12.00 

12.00 

12.00 

.00 

RATIO 

0.33 

10. A7 

10.09 

9.09 

259. 92 

000.00 

RANK 

» 

3 

2 

1 

4 

5 

AG 

2.00 

.05 

1.70 

1.40 

12.34 

.67 

sia 

12.00 

12.00 

12.00 

12.00 

12.00 

A. 00 

RATIO 

.17 

.07 

.14 

.12 

1.03 

.11 

RANK 

1 

4 

3 

5 

2 

AL 

2AO.OO 

317.00 

300.21 

297.57 

429.29 

.00 

sia 

12.00 

12.00 

12.00 

12.00 

12.00 

.00 

RATIO 

22.33 

2A.49 

a.  02 

24.00 

35.77 

000.00 

RANK 

» 

3 

2 

1 

4 

5 

CR 

1A.00 

17.29 

17. 1A 

17.10 

39.32 

.00 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

.00 

RATIO 

1.33 

1.44 

1.43 

1.42 

3.20 

000.00 

RANK 

* 

3 

2 

1 

4 

5 

CU 

10A.0O 

102.05 

100.  OA 

99.74 

100.09 

107.97 

sia 

12.00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

0.03 

0.50 

0.34 

0.31 

0.34 

9.00 

RANK 

4 

2 

1 

3 

5 

MG 

91.00 

170.50 

154.21 

140.12 

179.10 

.00 

sia 

12.00 

12.00 

12.00 

12.00 

12.00 

.00 

RATIO 

7.50 

14.21 

13.02 

12.34 

14.92 

000.00 

RANK 

II 

3 

2 

1 

4 

5 

NZ 

30.00 

34.07 

26.02 

31.70 

192.41 

.00 

sia 

12.00 

12.00 

12.00 

12.00 

12.00 

.00 

RATIO 

3.17 

2.91 

2.17 

2.64 

16.03 

000.00 

RAM 

3 

1 

2 

4 

5 

SI 

1A2.00 

149.09 

141.29 

166.29 

171.09 

103.04 

sia 

12.00 

12.00 

12.00 

12.00 

12.00 

6.00 

RATIO 

13.50 

12.42 

11.77 

13.06 

14.32 

17.17 

RAM 

2 

1 

3 

4 

5 

TI 

135.00 

142.77 

364.12 

141.42 

134.41 

141.47 
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SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

11.25 

11.90 

30.34 

11.  T9 

11.20 

11.79 

RANK 

4 

5 

2 

1 

3 

DATAAL 

PPM404 

FE 

30.00 

38.60 

92.84 

73.23 

1614.32 

.00 

SIZE 

6.00 

6.00 

6.00 

6.00 

6.00 

.00 

RATIO 

5.00 

6.43 

15.47 

12.21 

269.05 

888.00 

RANK 

* 

1 

3 

2 

4 

5 

A6 

9.00 

7.92 

6.83 

8.49 

98.37 

.00 

SIZE 

6.00 

6  •  Cv 

6.00 

6.00 

6.00 

.00 

RATIO 

1.50 

1.32 

1.14 

1.42 

16.39 

888.00 

RANK 

2 

1 

3 

4 

5 

AL 

231.00 

427.14 

240.35 

327.70 

966.05 

.00 

SIZE 

6.00 

6.00 

6.00 

6.00 

6.00 

.00 

RATIO 

38.50 

71.19 

40.06 

54.62 

161.01 

888.00 

RANK 

* 

3 

1 

2 

4 

5 

CR 

1.00 

.97 

2.22 

1.10 

10.77 

.00 

SIZE 

6.00 

6.00 

6.00 

6.00 

6.00 

.00 

RATIO 

.17 

.16 

.37 

.18 

1.80 

888.00 

RAfK 

1 

3 

2 

4 

5 

CU 

113.00 

116.47 

125.29 

127.72 

238.24 

.00 

SIZE 

6.00 

6.00 

6.00 

6.00 

6.00 

.00 

RATIO 

18.83 

19.41 

20.88 

21.29 

39.71 

888.00 

RANK 

* 

1 

2 

3 

4 

5 

MG 

47.00 

44.59 

45.87 

46.35 

46.32 

.00 

SIZE 

6.00 

6.00 

6.00 

6.00 

6.00 

.00 

RATIO 

7.83 

7.43 

7.64 

7.73 

7.72 

888.00 

RAT* 

1 

2 

4 

3 

5 

NI 

11.00 

8.71 

5.84 

6.41 

145.33 

.00 

SIZE 

6.00 

6.00 

6.00 

6.00 

6.00 

.00 

RATIO 

1.83 

1.45 

.97 

1.07 

24.22 

888.00 

RAMC 

3 

1 

2 

4 

5 

SI 

36.00 

53.14 

42.51 

69.86 

820.92 

.00 

SIZE 

6.00 

6.00 

6.00 

6.00 

6.00 

.00 

RATIO 

6.00 

8.86 

7.08 

11.64 

136.82 

888.00 

RANK 

» 

2 

1 

3 

4 

5 

TI 

32.00 

26.54 

25.41 

26.86 

25.17 

.00 

SIZE 

6.00 

6.00 

6.00 

6.00 

6.00 

.00 

RATIO 

5.33 

4.42 

4.23 

4.48 

4.20 

888.00 

RAMC 

3 

2 

4 

1 

5 

DATAAL  PPM40S 


FE 

239.00 

223.84 

175.23 

193.22 

4758.95 

.00 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

.00 

RATIO 

19.92 

18.65 

14.60 

16.10 

396.58 

888.00 

RAMC 

3 

1 

2 

4 

5 

AG 

24.00 

19.51 

18.68 

19.07 

23.69 

14.31 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

2.00 

1.63 

1.56 

1.59 

1.97 

1.19 

RAMI 

4 

2 

3 

5 

1 

AL 

120.00 

136.21 

133.36 

135.06 

138.79 

133.60 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

10.00 

11.35 

11.11 

11.25 

11.57 

11.13 

RANK 

M 

4 

1 

3 

5 

2 
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at 

3.00 

4.47 

4.41 

4.10 

14.73 

3.51 

size 

IE. 00 

12.00 

12.00 

12.00 

12.00 

4.00 

RATIO 

.23 

.50 

.41 

.51 

1.23 

.54 

RANK 

ft 

3 

1 

2 

5 

4 

CO 

47.00 

41.02 

47.51 

44.43 

40.57 

44.51 

size 

12.00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

3.42 

3.42 

3.44 

3.07 

4.05 

4.13 

RANK 

1 

3 

2 

4 

5 

MB 

244.00 

254.14 

220.37 

144.30 

213.04 

74.34 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

4.00 

RATIO 

22.17 

21.10 

10.34 

14.14 

17.02 

12.34 

RAW 

5 

4 

2 

3 

1 

MI 

24.00 

32.24 

0.74 

20.00 

147.12 

.00 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

.00 

RATIO 

2.17 

2.44 

.73 

1.47 

13.43 

000.00 

RANK 

3 

1 

2 

4 

5 

SI 

45.00 

71.11 

71.03 

71.21 

05.72 

40.57 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

5.42 

5.43 

5.42 

5.43 

7.14 

5.05 

RAW 

3 

2 

4 

5 

1 

TI 

40.00 

44.50 

45.33 

43.02 

44.24 

25.10 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

4.00 

RATIO 

3.33 

3.71 

3.70 

3.45 

3.44 

4.20 

RANK 

ft 

3 

4 

1 

2 

5 

DATAAL 

PPM400 

FE 

444.00 

547.44 

343.14 

415.10 

17220.25 

.00 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

.00 

RATIO 

34.00 

45.42 

30.24 

34.40 

1435.02 

888.00 

RANK 

3 

1 

2 

5 

4 

AO 

33.00 

20.00 

24.34 

24.41 

37.04 

.00 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

.00 

RATIO 

2.75 

2.34 

2.14 

2.24 

3.04 

888.00 

RAW 

3 

1 

2 

4 

5 

Al 

74.00 

43.13 

50.35 

41.43 

144.71 

10.28 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

4.00 

RATIO 

4.17 

5.24 

4.04 

5.12 

12.48 

1.71 

RAW 

4 

2 

3 

5 

1 

CR 

2.00 

3.40 

3.37 

3.34 

44.84 

.00 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

.00 

RATIO 

.17 

.24 

.20 

.20 

8.07 

888.00 

RANK 

ft 

3 

2 

1 

4 

5 

CU 

32.00 

35.41 

41.01 

42.44 

43.21 

45.57 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

2.47 

2.45 

3.42 

3.58 

3.40 

3.80 

RANK 

ft 

1 

2 

3 

4 

5 

MB 

41.00 

03.41 

45.54 

44.42 

280.57 

.00 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

.00 

RATIO 

7.50 

4.45 

7.47 

7.87 

23.38 

888.00 

RAW 

1 

3 

2 

4 

5 

HI 

30.00 

20.04 

0.25 

14.48 

478.03 

.00 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

.00 

RATIO 

2.50 

2.40 

.44 

1.41 

34.84 

888.00 

RANK 

3 

1 

2 

4 

5 
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sx 

132.00 

329.10 

306.22 

305.16 

316.09 

.00 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

.00 

RATIO 

11.00 

27.92 

25.52 

25.43 

26.34 

880.00 

RANK 

« 

4 

2 

1 

3 

5 

TI 

35.00 

32.77 

27.75 

29.39 

27.90 

28.11 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

2.92 

2.73 

2.31 

2.45 

2.33 

2.34 

RANK 

5 

1 

4 

2 

3 

OATAAL 

PPM409 

FE 

915.00 

896.53 

673.04 

788.52 

2643.86 

.00 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

.00 

RATIO 

76.25 

74.71 

56.09 

65.71 

220.32 

888.00 

RANK 

3 

1 

2 

4 

5 

AO 

41.00 

147.67 

5523.28 

292.54 

90.84 

7.25 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

6.00 

RATIO 

3.42 

12.31 

460.27 

24.38 

7.57 

1.21 

RAFK 

3 

5 

4 

2 

1 

AL 

277.00 

283.60 

265.18 

264.74 

255.95 

265.68 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

23.08 

23.63 

22.10 

22.06 

21.33 

22.14 

RANK 

5 

3 

2 

1 

4 

CR 

3.00 

4.47 

5.43 

5.10 

26.16 

.00 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

.00 

RATIO 

.25 

.37 

.45 

.42 

2.18 

888.00 

RANK 

If 

1 

3 

2 

4 

5 

CU 

134.00 

146.62 

157.18 

154.99 

218.25 

153.89 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

11.17 

12.22 

13.10 

12.92 

18.19 

12.82 

RAML 

If 

1 

4 

3 

5 

2 

MG 

73.00 

77.80 

78.61 

76.53 

95.99 

66.93 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

6.00 

RATIO 

6.08 

6.48 

6.55 

6.38 

8.00 

11.15 

RANK 

If 

2 

3 

1 

4 

5 

NI 

3.00 

2.85 

73.03 

6.58 

97.48 

7.02 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

6.00 

RATIO 

.25 

.24 

6.09 

.55 

8.12 

1.17 

RANK 

1 

4 

2 

5 

3 

SI 

525.00 

.00 

216.03 

.00 

5087.06 

.00 

SIZE 

12.00 

.00 

12.00 

.00 

12.00 

.00 

RATIO 

43.75 

888.00 

18.00 

888.00 

423.92 

888.00 

RANK 

3 

1 

4 

2 

5 

TI 

29.00 

25.84 

26.30 

23.93 

26.04 

21.23 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

6.00 

RATIO 

2.42 

2.15 

2.19 

1.99 

2.17 

3.54 

RANK 

2 

4 

1 

3 

5 

OATAAL 

PRM412 

FE 

41.00 

50.60 

142.39 

76.18 

1098.40 

85.99 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

6.00 

RATIO 

3.42 

4.22 

11.87 

6.35 

91.53 

14.33 

RAFK 

* 

1 

3 

2 

5 

4 

AO 

6.00 

5.50 

3.93 

4.56 

11.75 

1.68 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

.50 

.46 

.33 

.38 

.98 

.14 
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RAM 

4 

2 

5 

5 

1 

AL 

65.00 

241.77 

557.54 

260.60 

55.54 

.00 

SIZE 

It.  00 

It. 00 

12.00 

12.00 

6.00 

.00 

RATIO 

7.06 

20. 15 

29.78 

25.40 

8.89 

666.00 

RANK 

ft 

2 

4 

5 

1 

5 

CR 

5.00 

2.55 

1.94 

2.47 

15.26 

.25 

SIZE 

It.  00 

12.00 

12.00 

12.00 

12.00 

6.00 

RATIO 

.ts 

.21 

.16 

.21 

1.27 

.04 

RANK 

4 

2 

5 

5 

1 

CU 

55.00 

55.62 

47.50 

44.54 

51.29 

47.56 

SIZE 

It. 00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

2.75 

2.44 

5.94 

5.70 

4.27 

5.95 

RANK 

» 

1 

5 

2 

5 

4 

HO 

55.00 

42.21 

45.64 

44.11 

165.60 

16.44 

SIZE 

It.  00 

12.00 

12.00 

12.00 

12.00 

6.00 

RATIO 

4.42 

5.52 

5.64 

5.66 

15.65 

5.07 

RAM 

2 

5 

4 

5 

1 

NI 

25.00 

14.25 

609.15 

19.54 

96.26 

22.57 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

6.00 

RATIO 

1.92 

1.55 

67.45 

1.61 

6.02 

5.76 

RANK 

1 

5 

2 

4 

5 

SI 

680.00 

1541.92 

.00 

.00 

5279.27 

.00 

SIZE 

12.00 

12.00 

.00 

.00 

12.00 

.00 

RATIO 

75.55 

128.49 

666.00 

688.00 

275.27 

888.00 

RANK 

ft 

1 

5 

4 

2 

5 

TI 

21.00 

21.51 

20.92 

20.66 

21.59 

20.88 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

1.75 

1.78 

1.74 

1.72 

1.78 

1.74 

RANK 

4 

5 

1 

5 

2 

DATAAL  PPH413 

fl 

320.00 

298.66 

2303.91 

231.51 

966.67 

147.91 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

26.67 

24.89 

191.99 

19.29 

62.22 

12.33 

RANK 

3 

5 

2 

4 

1 

AG 

17.00 

17.38 

15.35 

16.32 

21.01 

14.15 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

1.42 

1.45 

1.28 

1.36 

1.75 

1.18 

RANK 

4 

2 

3 

5 

1 

AL 

124.00 

267.50 

136.11 

14.47 

16.94 

16.80 

SIZE 

12.00 

12.00 

12.00 

6.00 

6.00 

6.00 

RATIO 

10.33 

22.27 

11.51 

2.41 

2.62 

2.80 

RANK 

5 

4 

1 

3 

2 

CR 

7.00 

4.46 

20.26 

4.47 

15.25 

1.18 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

6.00 

RATIO 

.58 

.57 

1.69 

.37 

1.27 

.20 

RANK 

5 

5 

2 

4 

1 

CU 

81.00 

65.56 

66.28 

61.92 

99.51 

57.06 

SIZE 

12.00 

12.00 

12.00 

12.00 

12.00 

12.00 

RATIO 

6.75 

5.45 

5.69 
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